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ABSTRACT 

RVSAO is a set of programs to obtain redshifts and radial velocities from digital 
spectra. RVSAO operates in the IRAF(Tody 1986, 1993) environment. The heart of 
the system is xcsao, which implements the cross-correlation method, and is a direct 
descendant of the system built by Tonry and Davis (1979). emsao uses intelligent 
heuristics to search for emission lines in spectra, then fits them to obtain a redshift. 
sumspec shifts and sums spectra to build templates for cross-correlation, linespec 
builds synthetic spectra given a list of spectral lines, bcvcorr corrects velocities for the 
motion of the earth. We discuss in detail the parameters necessary to run xcsao and 
emsao properly. 

We discuss the reliability and error associated with xcsao derived redshifts. We 
develop an internal error estimator, and we show how large, stable surveys can be used 
to develop more accurate error estimators. 

We develop a new methodology for building spectral templates for galaxy redshifts, 
using the new templates for the FAST spectrograph (Fabricant, et al, 1998) as an ex- 
ample. We show how to obtain correlation velocities using emission line templates. 
Emission line correlations are substantially more efficient than the previous standard 
technique, automated emission line fitting. 

Using this machinery the blunder rate for redshift measurements can be kept near 
zero; the automation rate for FAST spectra is ~95%. 

We use emsao to measure the instrumental zero point offset and instrumental sta- 
bility of the Z-Machine and FAST spectrographs. 

We compare the use of RVSAO with new methods, which use Singular Value Decom- 
position and x 2 fitting techniques, and conclude that the methods we use are either equal 
or superior. We show that a two-dimensional spectral classification of galaxy spectra can 
be developed using our emission and absorption line templates as physically orthogonal 
basis vectors. 



Subject headings: methods: data analysis 
techniques: radial velocities 
instrumentation: spectrographs 
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1. Introduction 

Radial velocities are, along with position and bright- 
ness, among the fundamental measured values of as- 
tronomy. Recent technical advances are substantially 
increasing our ability to aquire radial velocity data; 
in the decade of the 1990's the rate at which radial 
velocity measurements are taken will increase by two 
or three orders of magnitude. Substantial effort is re- 
quired for these data to be reduced and analyzed in 
an accurate and timely fashion; here we describe the 
current reduction methods which we have developed 
for use by the Center for Astrophysics radial velocity 
and redshift programs, as well as by others. 

Doppler (1841) understood that radial velocities 
would affect the color of stars (by analogy with the 
pitch of sound); Fizeau (1848, 1870) first recognized 
that this would mean a shift in the position of the 
Fraunhofer lines. Huggins (1868) made the first (vi- 
sual) attempt (in 1862) to observe the shifts. Vogel 
(1892) made the first accurate photographic measure- 
ments, and established most of the procedures neces- 
sary to calibrate and reduce the measurements of line 
positions to radial velocities. 

Correlation methods for obtaining radial velocities 
were first suggested by Fellgett (1953), who was influ- 
enced by radar studies during World War II. Griffin 
(1967) was the first to implement these techniques. 
Note that Griffin credits Evershed (1913) with in- 
venting the basic technique; Griffin further notes that 
Babcock (1955) had already built a similar instru- 
ment. Griffin's instrument performed analog corre- 
lations by physically shifting a template spectrum in 
the focal plane of the spectrograph, a technique which 
is still in heavy use today (e.g. Baranne, et al, 1979). 

Digital power spectrum techniques for the estima- 
tion of lag have long been known (e.g. Blackman and 
Tukey, 1958, and references therein). Their use first 
became practical with the advent of digital detectors, 
fast digital computers, and the FFT algorithm (Coo- 
ley and Tukey, 1965). 

Simkin (1974) first showed how Fourier techniques 
could be used to obtain radial velocities and veloc- 
ity dispersions from digital spectra. Several groups 
used power spectrum techniques to obtain velocity 
dispersions (see Sargent, et al 1977, and references 
therein) and obtained velocities as a byproduct of 
their analysis, but apparently the first use of digi- 
tal cross-correlation specifically to obtain radial ve- 
locities was by Lacy (1977) who did not use Fourier 



techniques, but used direct convolution with a digital 
mask, emulating Griffin's (1967) analog technique. 

Tonry and Davis (1979; hereafter TD79) studied 
the use of power spectrum techniques to obtain red- 
shifts from digital spectra, and demonstrated the ef- 
fectiveness of the method. TD79 invented the r statis- 
tic, which can be calibrated to give both the confi- 
dence and error of a measurement. The techniques 
and software described here are directly descended 
from the TD79 system; in September 1990 radial ve- 
locity reductions at the CfA were moved from the old 
Data General Nova computer where the TD79 system 
resided onto a Unix workstation and the IRAF (Tody, 
1986, 1993) environment. 

We began with the IRAF task XCOR by G. Kriss 
and routines from TD79, translated into Fortran by 
J. Tonry ( Tonry and Wyatt 1988| ); these were ex- 
tensively modified, and resulted in xcsao version 1.0 
(Kurtz, et al, 1992; Paper 1). Additionally we used 
algorithms from th e REDUCE/INTERACT system 
( Maker et al., 1982Q , as modified by J. Thorstensen. 
The emission line finding programs had a somewhat 
different history. When the radial velocity reduc- 
tions were moved from the Novas onto Unix the emis- 
sion line programs were implemented as stand-alone 
C programs translated from the FORTH of TD79 by 
W. Wyatt; work began in 1991 on a new IRAF task, 
resulting in emsao (Mink and Wyatt, 1995). 

The software described here has been used ex- 
tensively; examples include the redshift surveys of 
Huchra, et al (1995); da Costa, et al (1994), Sheet- 
man, et al (LCRS:1996), Vettolani, et al (1997), and 
Geller, et al (1997). Stellar use revolves around the 
CfA Digital Speedometry program (Latham, 1985) 
and its many projects and collaborations. Latham 
(1992) lists several of these projects. 

Nordstrom, et al (1994) have described the tech- 
niques used by the CfA stellar group in detail; we 
will concentrate on issues related to galaxy redshift s, 
although some stellar data will be used in section || 
on error. We will conform to the convention that user 
settable parameters will be in italics, PARAMETERS 
in the SPECTRUM HEADER will be in CAPITALS, 
and IRAF tasknames will be in lowercase bold. 
Extensive on-line documentation, help files, and ex- 
amples, as well as the source code and execu tables 
can be found at http:/ /tdc-www. harvard.edu/ . 
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2. Practical use of xcsao 

xcsao is the heart of the RVSAO system. In this 
section we give an overview of the correlation system. 
We point out critical features which investigators need 
to consider when setting up the reductions for a new 
project, to maximize the efficiency of measurement, 
and minimize systematic errors. The RVSAO system 
consists of several IRAF tasks, including xcsao, the 
algorithmic details of each of them are described in 
appendix |a|. 

2.1. spectrum preparation 

Tokarz and Roll (1997) discuss the steps we take 
to obtain 1-D wavelength calibrated spectra suitable 
for redshift measurements. Once the 1-D spectra are 
in hand it is necessary to tell xcsao the wavelength 
range over which the data are good, using the param- 
eters stJambda and endJambda. Depending on the 
details of the instrumental set up substantial addi- 
tional error can occur if these parameters are not set, 
or are set incorrectly. For example, if a substantial 
portion of the spectrum is from a region of the de- 
tector with little or no sensitivity, or is from a region 
of the spectrum with poor sky subtraction and strong 
night sky lines, the final redshift will be compromised. 

2.2. continuum and emission line suppres- 
sion 

For all spectra the con tinuum must be removed; 
contpars (see section AJ3) performs these tasks here, 
and was only slightly modified from the IRAF ONED- 
SPEC conti nuum package, as implemented in the 
RV package flFitzpatrick 1993|) . 

The continuum can be subtracted or divided out. 
Subtraction preserves the correct relative amplitudes 
for the lines in the data, and thus the correct signal to 
noise behavior in the cross-correlation. Division pre- 
serves the correct equivalent widths of the lines, and 
thus the correct parameterization of the spectrum. 
For spectral classification studies division is preferred 
(Kurtz, 1982); for radial velocity measurements sub- 
traction is superior. Division by the continuum re- 
sults in amplified noise in the blue part (the low S/N 
part) of the spectrum. For moderate S/N spectra, 
typical of FAST, the difference between the two tech- 
niques is small, but subtraction shows smaller redshift 
residuals by about 25% ; while for very high S/N spec- 
tra on FAST, typical of our calibration spectra, the 
division method gives slightly smaller residuals. 



In normal use of xcsao only continuum subtraction 
is allowed, but if the parameter DIVCONT is set true 
(T) in the template spectrum header the continua of 
both the object and template spectra will be divided 
out. 

Many galaxy spectra show both emission and ab- 
sorption lines; in general the redshifts derived from 
the emission lines will be different from the absorp- 
tion line velocities. To obtain a correlation velocity 
from an absorption line template for a spectrum with 
strong emission lines it is necessary to suppress the 
emission lines; figure |l| shows a typical spectrum and 
its correlation function with one of our standard ab- 
sorption templates suppressing the emission lines and 
without suppressing them, the bottom panel shows 
the spectrum, smoothed and with the lines marked. 
The reduction with emission line suppression pro- 
duces a believable redshift, with an r (TD79) value 
of 4.50, the reduction without it, where the r value is 
1.88, does not. 

Removing emission lines before correlating with an 
absorption line template was a routine feature of the 
TD79 software; in xcsao we have extended and gen- 
eralized the procedure. Both emission and absorption 
lines may now be removed, and the process may be 
controlled by keywords in the template header. Be- 
cause we now obtain emis sion line velocities using cor- 
relation methods (section |4.l|) , the same emission line 
suppression cannot be used for every template. 

For routine redshift reductions at the CfA we use 
a subset of the xcsao capabilities, with emission lines 
replaced by the continuum when we correlate against 
an absorption line template, and the absorption lines 
replaced by the continuum when we correlate against 
an emission line template. The exact parameters used 
are stored in each template, normally two sigma vari- 
ations above (below) the continuum are sufficient to 
remove emission(absorption) lines, with the number 
of iterations and growing parameter as set in cont- 
pars (sect ion ^5|). 

The emission(absorption) line suppression is con- 
trolled by the parameters S-emchop (t-emchop) for the 
object (template) spectra. Other parameters are used 
as well, and when control is given to the template or 
object spectrum header the interactions can be com- 
plex. They are more fully described in section |A.l] , 

In addition xcsao permits the user to replace spe- 
cific regions of the spectrum with a simple linear ap- 
proximation to the continuum. This feature is typi- 
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Fig. 1. — The effect of emission line suppression on 
absorption line correlation redshifts. The upper panel 
shows the observed spectrum; the next panel down 
shows the correlation function with an absorption line 
template after emission line clipping; the next down 
shows the correlation function when the emission lines 
are not clipped; the bottom panel shows the observed 
spectrum after smoothing, and with the main lines 
marked. 

cally used to prevent poor subtraction of the bright 
night sky lines from compromising the results. 

2.3. apodization, zero padding, Fourier fil- 
tering 

Once the continuum has been removed the spec- 
tra are apodized, zero padded, and bandpass filtered. 
Each of these operations has a substantial effect on 
the final results. 

Apodization is the simplest; essentially the goal is 
to remove any ringing in the Fourier transform, by 
forcing the ends of the spectrum smoothly to zero, 
while suppressing as little of the actual data as possi- 
ble. Our apodization is performed by a cosine taper 
function, which begins a symmetric set percentage of 
the spectrum from the ends. For both FAST spec- 
tra and the earlier Z-Machine (Latham, 1982) spectra 
0.05 is a reasonable value, i.e. the taper begins 5% 



from the ends (section A.l). 

Zero padding the spectrum is intended to remove 
any artifacts caused by computation of the correla- 
tions in Fourier space, thus using a circular convolu- 
tion. The primary artifact which the zero padding 
removes is the confusion of Ha with OII3727 due 
to the wrap around of the convolution. The zero 
padding has two side effects which must be consid- 
ered. First, the relation of the TD79 r statistic with 
error is changed, although the actual calculated error 
remains correct (section ||). Second, the envelope of 
noise fluctuations of the correlation function, which 
is flat in the non-zero padded case, is, in the zero 
padded case, a symmetric linear function of the num- 
ber of overlapping non-zero pixels, and is maximum 
at the redshift of the template spectrum. For the case 
of low S/N spectra where the redshift is substantially 
different from the template's redshift this structure 
in the noise can result in the wrong correlation peak 
being chosen. Zero padding can be controlled via the 
parameter file, or the template header. We only zero 
pad spectra when correlating against the emission line 
template. 

The design of the Fourier bandpass filter is criti- 
cal to the optimal measurement of redshifts. As with 
the apodization we use a cosine taper to suppress the 
ends of the (in this case) Fourier spectrum. Several 
other filter techniques were tried (e.g. Oppenheim 
and Schafer, 1975) but no difference was seen for any 
reasonable choice of taper function (a sharp cutoff is 
not reasonable because of Gibbs ringing). In addi- 
tion we tested a spectral weighting function shown 
by Hassab and Boucher (1979) to produce the maxi- 
mum likelihood estimator for the lag (radial velocity) 
in the limit of infinitely wide spectra; this weighting 
function had no positive effect on our results, and we 
have not implemented it. 

Removing high spatial frequency information, via 
a high- stop Fourier filter, is intended to increase the 
S/N by removing information which contains more 
noise than signal. The design question is where to 
set the high frequency turnoff; the method we use 
is to examine sets of high S/N calibration spectra, 
e.g. our nightly exposures of NGC4486b. We cor- 
relate each of these against the best match template 
(in this case NGC7331), using a high pass Fourier 
filter which filters out ALL the low frequency infor- 
mation, leaving only the high frequency noise. If the 
turn-on frequency is set too high, the correlations all 
give an incorrect redshift; if the turn-on frequency is 
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set too low, all correlations give the correct redshift. 
We choose the turn-on frequency where half the red- 
shifts are correct and half incorrect, and set this turn- 
on-frequency to the turn-off frequency for our high- 
stop filter. This procedure gives a turn-off frequency 
approximately equal to that obtained by the "opti- 
mal filter" method of Brault and White (1971) which 
chooses that point where the power of the signal is 
twice that of the photon noise. This frequency is less 
than half that which corresponds to the projected slit 
width of the FAST. 

The high-stop filter is only used for absorption line 
spectra. For emission line spectra the redshifts are 
seriously degraded if the high spatial frequencies are 
removed from the data. The key word FI-FLAG in 
the template header controls the implementation of 
the high-stop filter. In addition, as it is possible to 
pre- filter the templates, to save unnecessary comput- 
ing, this flag also con trols whether and how to filter 
the templates. Figure A31 shows all the possibilities. 

Removing low spatial frequency information, via 
a low-stop Fourier filter, is intended to remove any 
residual large-scale systematics which remain follow- 
ing the continuum suppression. Essentially this can 
be viewed as a second continuum removal, equiva- 
lent to the continuum removal technique of LaSala 
and Kurtz (1985), thus giving what Kurtz and LaSala 
(1991) call a "reflattened" spectrum. 

The design question is where to put the low fre- 
quency turn-on point. The difficulty with making this 
decision is that there is no point where excluding all 
information with higher spatial frequencies does not 
result in a redshift (i.e. even the lowest spatial fre- 
quencies still contain accurate redshift information), 
and there is no reasonable point where suppressing 
more low frequency information does not result in 
lower residuals for high S/N sets of spectra, such as 
our set of NGC4486b spectra. 

We therefore set the low frequency turn-on point 
by a simple heuristic. We estimate the scale in wave- 
length of the broadest spectral feature useful in esti- 
mating a redshift, in the case of FAST galaxy spec- 
tra this is the change in the slope of the contin- 
uum around the Call H+K lines, and we suppress 
those spatial frequencies which correspond to twice 
this scale, or greater. 

The remaining design decision for the Fourier fil- 
ter is the width of the turn-on and turn-off ramps. 
It may be expected that this is only important at 



the low frequency turn-on point, as the power there 
is typically two orders of magnitude above the high 
frequency turn-off point. The problem is that if the 
turn-on is too sharp, Gibbs ringing will be introduced 
into the data. A full turn-on width of 1.5% of the 
width of the power spectrum is sufficient to amelio- 
rate this effect. 

The exact filter implemented, especially the exact 
implementation of the low-stop filter, affects the re- 
sulting redshifts, their errors, and the relation of the 
TD79 r statistic with their errors. For example, for 
the set of NGC4486b spectra, the mean redshift ob- 
tained using only the lowest spatial frequencies which 
we include in our standard filter differs from the mean 
redshift obtained by only using the highest included 
spatial frequencies by 61/cm/s, and different reason- 
able choices for the Fourier filter can give redshifts 
which differ in the mean by 10km/ s. These differ- 
ences may be compared with a typical variation about 
the mean of 15km/ s(la). The sign and amplitude 
of this effect changes with each object- template pair, 
NGC4486b vs. NGC7331 is a typical result. In ad- 
dition the relation of the TD79 r statistic with error 
depends on the filter (section ||) . 

2.4. cross-correlation, rebinning, and red- 
shift evaluation 

The cross correlation is the normal product of the 
Fourier transform of the object spectrum with the 
conjugate of the transform of the template spectrum, 
as described in TD79. 

The object spectrum and the template spectra 
need to be pairwise rebinned to have a common dis- 
persion. The number of bins nbins is set by the user 
and must be a power of two; we recommend that nbins 
always be larger than the number of observed pixels. 
The spectral region rebinned is set to obtain the max- 
imum overlap between the template spectrum and the 
portion of the object spectrum between stJambda and 
enddambda in the rest frame. On the first pass the 
rest frame is determined by a user guess to the red- 
shift (czguess), or from a previous reduction (section 
|A.ip . On subsequent passes (if nzpass>0) the rest 
frame is determined from the redshift obtained in the 
previous pass. We recommend that czguess be set to 
the approximate redshift expected (normally a better 
guess than zero), and nzpass=2. 

Next the correlation peak is determined and fit 
(section [A. if) . The type of fit (pkmode) has little effect 
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on the result, but the amount of the peak which is fit, 
pkfrac, is critical. The fit is performed from the top of 
the peak down to where the peak is pkfrac of the max- 
imum, thus more of the peak is fit if pkfrac=0.5 than 
if pkfrac— 0.7. Because of side-lobes in the correlation 
function due to the proximity of Nil to Ha, emission 
line correlation peaks cannot be fit as far down the 
peak as absorption line correlations. The template 
header parameter PEAKFRAC overrides pkfrac on a 
template by template basis; we use this to set the fit 
parameters for the emission line templates. 

3. Error Analysis 

There are three main questions concerning the out- 
put of xcsao: are the results reliable, can the process 
be automated? what is the size and nature (random 
or systematic) of the error? and, does what is mea- 
sured correspond to the physical property the inves- 
tigator wants to measure? 

3.1. reliability 

To answer the first two of these questions we use 
a dataset designed for this purpose; it contains 626 
pairs of spectra observed with the FAST spectrograph 
between 1994 and 1996. Each pair consists of two 
independent observations of the same object; about 
half of these were observed to calibrate the velocity 
errors for the 15R survey (Geller, et al. 1998), and 
about half are spectra which were below the quality 
standard, and required a second integration (these 
would be summed in our normal reductions, but not 
here). All these spectra have been subjected to our 
normal processing (Tokarz and Roll, 1996), and thus 
have had most of the cosmic rays removed by a labor 
intensive process. 

As discussed by TD79, the r statistic can be cal- 
ibrated as a confidence measure. We prefer to cali- 
brate it empirically, rather than use the prescription 
in TD79. Figure |2| shows the results of correlating 
each of the 1252 spectra against each of two tem- 
plates. Plotted are the absolute value of the veloc- 
ity difference between the two observations versus the 
minimum r value for each pair. Each pair appears on 
the graph twice. The open circles are measurements 
using the NGC7331 template, and the filled triangles 
use an emission line template, emtemp. 

Figure [| shows two groupings, those with the 
absolute velocity difference Av ~ 300km/ s; which 
we will assume are reliable observations (differences 
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Fig. 2. — The blunder diagram, see text for discus- 
sion The thick dotted lines are at r — 3, and velocity 
differences of 300 km/s; the thin dashed line is at 
r = 2; both axes are on a log scale. 

~ 300km /s are consistent with the expected ran- 
dom errors), and those where the velocity difference 
is greater, which we will assume unreliable. For the 
moment we will adopt an r value of 3.0 (r m j n ), above 
which we expect the velocity determination to be re- 
liable. This is lower than is typically used for FAST 
reductions. 

We would then expect no points in the upper right 
quadrant of the plot, where r > r m ^ n , and the Av > 
300 km/s] however there are, 15 points in what we 
will call the "blunder" region. 

We will examine each point to determine which 
rules would catch the blunders in a fully automated 
reduction. The two nearly co-incident circles (cir- 
cles are reductions with the NGC7331 template) at 
r ~ 3.2 and velocity difference about 50000/cm/s are 
both objects with strong emission lines. It is known 
that the NGC7331 template systematically returns 
velocities of 49000/cm/ s for some emission line ob- 
jects, so along with the high r value emission line 
velocity these measures could be discarded automat- 
ically. The circle at r ~ 3.1, Av ~ 18000/cra/s is 
also an object with strong emission lines. A simple 
rule which requires that spectra with discordant, but 
otherwise valid velocities be checked manually would 
catch this, but the reduction could not be fully auto- 
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matic. 

The four triangles (triangles are the emission line 
template) between r ~ 3.2 and r ~ 4.2 having a 
Av ~ 948/cm/s all have good absorption line veloci- 
ties, and would be caught by the rule of discordance. 
They would be caught by another rule, however, one 
which does not require that the absorption line veloc- 
ity be "good;" they are all spectra where Nil 6583A is 
stronger than Ha, and the difference with the absorp- 
tion line templates is about 948/cra/s. We adopt the 
rule that all emission line velocities which differ from 
an absorption line velocity (even if one with a low r 
value) by about 948/cm/s must be checked manually. 

The circle near r ~ 4.5 and Av ~ 350km/ s is prob- 
ably not a real blunder. Examination of the POSS 
prints shows that the object has two nucleii with 4" 
separation. We assume that the velocity difference is 
real, and that the two observations of this object each 
correspond to a different nucleus. 

All of the seven remaining objects in the "blunder" 
region of the plot are emission line velocities. One has 
the night sky line at 5577A mistaken for Oil 5007A; 
we can eliminate this error by either turning the bad- 
lines removal feature on to replace the region around 
5577A with the continuum, or by adopting the rule 
that all emission line redshifts near 34152/cra/s be ex- 
amined manually. 

The remaining six are all spectra contaminated by 
cosmic rays. Five of these would be tagged by the 
rule of discordance, and the sixth has an absorption 
line r value of 2.77, so it would be tagged by only a 
slightly more stringent rule of discordance; it cannot 
be assumed, however, that the cosmic ray problem 
can be solved by looking for absorption lines. We 
therefore require that emission line red shifts must all 
be checked using emsao (section |A.2|) , and that at 
least four lines must be found which correspond to 
the correlation velocity, and at least two must be fit; 
"blunders" are made when only three lines are found, 
or only one line is fit. Using that criterion all six 
spectra would be tagged for visual inspection as well 
as 30 of the 178 emission line velocities (r > 3) which 
do not have a confirming absorption line velocity with 
r > 3. 

For the 610 objects (of 626 total objects observed) 
where at least one of the two different template re- 
ductions gave a result with r > 3, xcsao yielded the 
correct result with no further problem for 595. Of the 
remaining 15 spectra 12 are easily discovered because 



two valid redshift measures disagree; one is probably 
caused by source confusion on the sky, and two are 
found by emsao. 

Looking at figure ^, it is clear that many spectra 
where 2 < r < 3 do indeed give the correct redshift. 
Of the 16 objects which have neither emission nor 
absorption reduction with r > 3 four have both emis- 
sion and absorption redshifts equal (within normal er- 
rors) and could be accepted (we do not currently do 
this) . Also 66 spectra where the emission line r value 
is > 3 have confirming absorption line velocities with 
2 < r < 3; if these are assumed correct (which we also 
do not currently do) , then the number of emission line 
spectra which must be visually inspected after emsao 
would drop from 30 to 19. This would bring the total 
number which require visual inspection to 31, or 5% 

With the aid of emsao for quality control, and par- 
tial manual reduction of 31 spectra, xcsao obtained 
the correct redshift for all 614 objects which yielded 
redshifts, save for the one object which was probably 
an observational error. 

A second experiment was made using 8606 emis- 
sion line spectra from the Z-Machine archive, which 
had r values with the emission line template above 
3. The results of the correlation with the emission 
line template were compared with the stored redshift 
in the archive, which was obtained by manually fit- 
ting the emission lines with a precursor program to 
emsao. After sifting the results using rules like those 
described above fifteen spectra (0.2% ) had the wrong 
redshift; essentially all these spectra were the victims 
of very poor sky subtraction. This may be compared 
with the twenty-four spectra where the redshifts were 
incorrectly listed in the archive. The Z-Machine spec- 
tra were substantially noisier than the FAST spectra; 
nearly 15% failed the sifting and would have had to 
be manually reduced. 

Used carefully the RVSAO suite provides redshifts 
with a very low blunder rate. The automation rate 
obtainable with RVSAO is strongly affected by the 
S/N of the observations. Absorption line objects must 
be observed long enough to have a fully reliable ab- 
sorption correlation velocity (we currently use r ~ 4, 
which is conservative) or a confirming weak emission 
velocity. Emission line spectra must have a confirm- 
ing weak absorption redshift, or be based on at least 
four lines, and at least two of the four must be fit by 
emsao (section flT|) . 
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3.2. error estimation 

Besides estimating the redshift of a spectrum xc- 
sao also estimates the error in the redshift. The 
error estimator can be derived analytically follow- 
ing the discussion in section IILc.i of TD79 with the 
additional assumption of sinusoidal noise, with the 
halfwidth of the sinusoid equal to the halfwidth of 
the correlation peak. The derived error estimator is: 



w 



error 



8(l + r) 

where error is the error in a single velocity mea- 
surement by xcsao, w is the FWHM of the correlation 
peak, and r is as defined in TD79. 

While the assumption of sinusoidal noise with 
halfwidth equal to the correlation peak's halfwidth is 
reasonable, there is no compelling argument for this 
assumption. Therefore it is necessary to demonstrate 
the effectiveness of the approximation by experiment. 

We will use four datasets to examine the behavior 
of the error estimator: the 610 duplicates described 
above; the 8606 Z-Machine emission line spectra de- 
scribed above; 7810 synthetic spectra, each identically 
Poisson sampled from a 45A section of a model atmo- 
sphere for a 5500K dwarf star (Kurucz, 1992), taken 
from the set of synthetic stellar templates used by 
the CfA digital speedometry program (Morse, et al 
1991; Nordstrom, et al 1994); 50000 synthetic spec- 
tra, using the same 5500K dwarf star template, each 
with a different number of simulated photons (we con- 
fine ourselves to using the 49880 spectra which, when 
correlated against the synthetic template, achieved 
r > 3). 

First we will look at the set of duplicate spectra. 
We will limit ourselves to cases where both reductions 
have r values > 3.5. This is 400 ±15 spectra for 
the two absorption line combinations, and 297 for the 
emission line comparison. 

Figure [| shows a typical result. The solid line 
shows a histogram of the absolute values of the differ- 
ences between two observations of the same object, 
both reduced in the same way using the NGC7331 
template, and divided by the sum in quadrature of 
the errors calculated by xcsao for the reductions. The 
dotted histogram is the expected Gaussian distribu- 
tion; it is clearly broader than the data, xcsao over- 
estimated the error by ~ 20% . 

Figure || is similar to Figure |[ Here the ztemp 
template was used on the FAST data; ztemp is a 
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Fig. 3. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the fn7331temp template, 
divided by the xcsao error estimate. The dotted his- 
togram is the expected Gaussian distribution. 



combination of bright galaxy spectra taken with the 
Z-Machine, and in use at the CfA since the days of 
TD79. ztemp has a restricted wavelength coverage 
(AA 4500 6200 A) compared with the FAST spectra, 
has a different resolution, and has different residual 
systematics. In Figure^ the (dotted) Gaussian is nar- 
rower than the (solid) data, xcsao underestimated 
the error by ~ 20% . 

Figure ^| shows a similar set of histograms for the 
emission line template, emtemp. emtemp is a syn- 
thetic spectrum m ade b efore the creation of the line- 
spec task (section |A.3|) to match FAST emission line 
galaxy spectra. Here the solid line which represents 
the data cannot be transformed to match the dotted 
line expected histogram by any multiplicative process 
(1.2 would be the best multiplicative factor); it would 
still have more power in the tail. Adding 15 km/s in 
quadrature helps remove power from the tail. 

TD79, while giving a procedure to calculate the er- 
ror, suggest that in practice the error be calculated 
by calibrating k in the equation error = k/{l + r) 
using external comparisons. Paper 1 reiterates this 
suggestion, noting that the measurement of w has er- 
ror, but when all reduction parameters remain fixed, 
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Fig. 4. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the ztemp template, di- 
vided by the xcsao error estimate. The dotted his- 
togram is the expected Gaussian. 
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Fig. 6. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the fn7331temp template, 
divided by the k/(l + r) error estimate. The dotted 
histogram is the expected Gaussian. 
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Fig. 5. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the emtemp template, di- 
vided by the xcsao error estimate. The dotted his- 
togram is the expected Gaussian. 



w is essentially constant. Figure II of Paper 1 shows 
the effect of changing one of the reduction parameters 
(the low frequency roll off of the Fourier filter) , which 
substantially changes the relation of (1 + r) to error, 
while w scales correctly so that 8 (i^ r ) still tracks the 
error. 

Given a set of duplicate observations the constant 
k can be determined by internal comparisons. The 
procedure is simply to vary k until the expected dif- 
ferences histogram matches the measured one. For 
the present case we obtain /catgc7331 = 315km/ s, 
kztemp = 285/cra/s, and k emtemv = 245km/ s. Fig- 
ures H |t], and H show the distributions compared with 
the expected Gaussians; as expected, in all cases the 
fit is better than with the unmodified xcsao errors. 
For large observing programs with stable reduction 
procedures, we recommend using calibrated h/(l + r) 
relations to estimate the error. 

The duplicate spectra, as is typical of redshift sur- 
vey data, do not show a very large range of S/N, or 
r value. This is due to the fact that one normally 
observes long enough to get a desired S/N, and no 
longer. If we continue to restrict the duplicate pairs 
to those where both spectra achieved r > 3.5 (still 
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Fig. 7. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the ztemp template, di- 
vided by the k/(l + r) error estimate. The dotted 
histogram is the expected Gaussian. 
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Fig. 8. — The solid histogram is the distribution of 
velocity differences for duplicate observations of the 
same object, reduced with the emtemp template, di- 
vided by the k/{\ + r) error estimate. The dotted 
histogram is the expected Gaussian. 




Fig. 9. — Residuals in fit to Gaussian error distribu- 
tion for pairs reduced with the emtemp template. The 
contours represent lines of equal residual, position in 
the x,y space represents admixtures of the statistical 
and systematic error models; see text. 

lower than is normally required for a FAST redshift) 
then it is not possible here to use the goodness of fit 
to a Gaussian to prove that the error calculated by 
using 1+r is any better than a constant error, for each 
template. If the error were dominated by systemat- 
ics one would expect the error to be approximately 
constant. 

We define a measure of error 



/ k 2 

/ r,2 _i_ statistical 

V ystematic (1 + ri) 2 



k 2 

statistical 

(l+r 2 ) 2 



We then vary the values of k 



systematic 



and k< 



tatistical 



and calculate the residuals of the fits to a Gaussian. 
Note that k sys t e matic is y/2 times the error in a sin- 
gle measurement. Figure [5] shows the results for the 
duplicate pairs reduced with emtemp. The contours 
represent lines of equal residuals in the fit to a Gaus- 



sian, in the (k 



statistical ? ^systematic 



) space. The outer 



contour is 20% larger than the inner contour. Clearly 
either a fully systematic or a fully statistical error is 
consistent with the data. 

Figure [l0| shows the same diagram for the N7331 
template. The inner contour level here represents an 
absolute error half that of emtemp, and the outer level 
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Fig. 10. — Residuals in fit to Gaussian error distri- 
bution for pairs reduced with the fn7331temp tem- 
plate. The contours represent lines of equal residual, 
position in the x,y space represents admixtures of the 
statistical and systematic error models; see text. 

is 33% larger than the inner level. Also here one can- 
not rule out either a fully systematic or a fully statis- 
tical error. 

In many cases it is not possible to reduce hun- 
dreds of duplicate measurements using exactly the 
same reduction procedures to obtain an improved er- 
ror estimator; in these cases the xcsao error estima- 
tor is a reasonable choice. The 20% systematic de- 
viations for the two absorption line galaxy templates 
are the largest we have seen, although Quintana, et 
al. (1996) suggest that for their data the error is un- 
derestimated by ~ 30%, by comparison with external 
measurements. While the xcsao error estimator dif- 
fers systematically from the true error for a particular 
combination of instrumental set-up, reduction proce- 
dure, and template we have not seen any trend for 
this to be a systematic over or under estimate. 

While the 8606 Z-Machine emission line spectra 
cannot be used to calibrate the error estimator, as 
to first order we are just comparing the differences 
in two different methods of fitting Ha in the same 
spectrum, we can use them to look at any differences 
in zero point, as a function of the fitting method. The 
mean difference between the correlation velocity and 
the velocity obtained by the semi-automated linefits 




Fig. 11. — The synthetic 5500K spectra. Above is 
the template spectrum(Kurucz, 1992), in the middle 
is a typical poisson sampled spectrum, and below is 
the same spectrum smoothed. 

is 0.56/cra/s ± .16. 

The 7810 synthetic spectra are each correlated 
against the template from which they were identically 
randomly Poisson sampled; thus there is no spec- 
tral type difference adding to the errors. Figure |ll 
shows the template, a typical sampled spectrum, and 
a smoothed version of the typical spectrum. By calcu- 
lating the RMS velocity about the expected velocity 
(0km/ s) we have a very accurate measure of the error 
in a single measurement; the mean error calculated by 
xcsao is 3% greater than this. The distribution of ve- 
locities is essentially Gaussian, and the deviation of 
the mean velocity is within 1 a of zero. 

Using the 7810 spectra we can ask the question: 
"how many independent measurements of a spectrum 
are required to obtain a better measure of the error in 
a single measurement than xcsao provides?" Nord- 
strom, et al (1994), on the basis of a study of echelle 
spectra of rotating F stars give this answer as ~ 7. 
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Fig. 12. — Variance in error estimation as a function 
of the number of measurements for the set of 7810 
synthetic spectra. The horizontal line is the variance 
(about the true value) of the xcsao error estimate. 
See text for details. 

Here we take as many independent sets of N spectra 
as exist in 7810 spectra, for each N we calculate the 
sample standard deviation about the sample mean, 
and we compare it with the known error gotten by 
using all 7810 velocities. The RMS of this difference 
is plotted as a function of N in figure [l^. Also plot- 
ted (as a straight line) is the RMS of the xcsao error 
estimate about the true error. In this ideal case more 
than 30 independent measures are required before a 
better error estimate is reached than the xcsao error. 

The 50,000 synthetic spectra, with 1 to 50,000 
counts, give similar results. The mean error is under- 
estimated here by 8% , and the mean velocity is 1.1 
a different from zero, using the 49880 spectra where 
r > 3. The efficiency of (1 +r) in estimating signal to 
noise is demonstrated in Figure [uj[ here we show the 
ratio of the square root of the number of counts to the 
best fit linear relation with (1 + r). The la scatter 
is ~ 12% , independent of N; this puts a limit on the 
inherent ability to estimate errors using (1 + r). 

3.3. systematics 

There are many factors which can cause the red- 
shifts and radial velocitys measured by xcsao to be 




Fig. 13. — The relation between r and y/N Each 
point represents the ratio of y/N with the best fit 
linear relation between y/N and 1 + r. 



other than those desired. Here we list several. 

1. Errors in the wavelength calibration. For ab- 
sorption line spectra, where the signal is averaged over 
all the lines in a complex spectrum, this effect should 
be tiny. For emission line spectra, where the signal 
comes from a couple of lines this effect could be as 
large as the error in the pixel to wavelength calibra- 
tion function. For typical FAST galaxy spectra this 
error is ~ 5/cm/s averaged over the entire wavelength 
range, it could be larger in small regions (section 4.3). 

2. Offsets in the calibration lamp illumination. 
The light from the calibration lamps does not follow 
the same exact optical path as the light from the sky. 
This can cause systematic errors in the wavelength 
scale (section ^|). 

3. Variatio ns w ith the Fourier filter. As noted 
above (section |2.3D different reasonable choices of the 
Fourier filter can change the measured velocity of a 
typical FAST galaxy spectrum by 10km/ s, compara- 
ble to the error in our highest S/N observations. Un- 
reasonable choices for the filter parameters can make 
more of a difference. 

4. Spectral type mismatch. It has long been known 
that there are systematic velocity effects when the 
template spectrum does not match the observed spec- 
trum. Nordstrom, et al. (1994), for example, show 
the effect of rotational velocity mismatch for their 
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echelle spectra. For typical FAST galaxy spect ra t his 
effect is ~ 
5. 



20km/ s, and is discussed in section [13 
HII regions in spiral galaxies. 



The rotation 

velocities of disks and the finite number of HII re- 
gions on the slit can combine to yield an Ha veloc- 
ity which is different from the mean velocity of the 
stars in the bulge; this is especially true for galaxies 
which are distorted. Thus, while the measurement er- 
ror in an emission line velocity for a particular galaxy 
may be substantially smaller than for the absorption 
line velocity, the systematic deviation from the de- 
sired quantity, the cosmological redshift, may be sub- 
stantially larger. 

6. Two-lined systems, xcsao assumes that the 
template spectrum is a reasonable spectral match to 
the observed spectrum. For the case of two (or more) 
lined spectroscopic binaries this condition is clearly 
violated. As demonstrated by Latham, et al. (1996) 
substantially improved results may be obtained by us- 
ing methods which explicitly model the two lined case, 
such as TODCOR (Zucker and Mazeh, 1994). xcsao 
may be used to obtain the input data for TODCOR. 

4. Templates 

Accurate redshifts require the existence of very 
high signal-to-noise templates, which have well de- 
termined velocities, and are good matches to the sci- 
entific program objects being measured. Here we de- 
scribe methods of creating and maintaining systems 
of templates. 

Templates must be very high signal to noise spec- 
tra; there are two basic ways to create them: 1) sum 
a number of observations; 2) build a computer model 
of a spectrum. The CfA Digital Speedometry group 
has, over the past decade, switched entirely from us- 
ing observed spectra to models (see Latham, et al, 
1996 and references therein for details). For galaxy 
redshift studies we use both techniques. 

The vast majority of nearby galaxies can be well 
matched by a typical absorption line spectrum (like 
NGC7331), by a typical emission line spectrum, or 
both. Unusual spectra, as one obtains for QSOs, RS 
strong galaxies, galaxies with extremely high or low 
internal velocity dispersions, etc. require special tem- 
plates; although in most cases special templates only 
lower the error in the redshift; the correct redshift is 
normally obtained using standard templates. 

Creating templates is an iterative process: good 
templates are required as a prerequisite for making 



better templates. New templates were made for FAST 
when it saw first light, and again when it received a 
new thinned CCD in September 1994. In 1997, suffi- 
cient new observations having been made, we were 
able to make a set of substantially improved tem- 
plates. 

4.1. emission line template 

Until 1995, redshifts for emission line galaxies were 
obtained with emsao and its precursor programs. 
Then a template was made by placing Gaussians with 
approximately correct line widths and line ratios at 
the emission line rest wavelengths. This template, 
emtemp, was tested against existing FAST observa- 
tions, and against the Z-Machine archive. Typical 
results are for the Z-Machine comparison with 8606 
emission line spectra reduced by hand: the la differ- 
ence is 13/cras -1 , about half of the calculated error for 
a typical Z-Machine emission line velocity. As noted 
above the zero point offset is 0.56km/ s ± .16, about 
a two hundredth of a pixel. The median difference 
of 3929 FAST spectra with r ernt emp > 3 and the H a 
velocity obtained by emsao is 0.24km/ s. 

Essentially for all spectra where emsao can obtain 
a redshift xcsao obtains a redshift using emtemp. For 
low S/N spectra xcsao plus emtemp is much more 
sensitive than emsao. For a sample of 2088 emis- 
sion line galaxy spectra taken with FAST only 42% of 
spectra where xcsao plus emtemp obtained a redshift 
with 3 < r emtemp < 4 could be reduced automatically 
using emsao, 64% for 4 < r ernte m P < 5, and 93% for 
r emtemp > 5. As noted above the systematic effect of 
cosmic rays places severe constraints on the unsuper- 
vised use of xcsao with an emission line template. 

To make a better emission line template we have 
taken 6498 FAST spectra where r ernte mp > 5, put 
them through emsao, and obtained 434 spectra where 
emsao found 9 or more lines. For these we measured 
the ratio of line heights with the SII 6731A line, and 
the line widths for each line. Then for each line we 
took the median of these quantities, and, along with 
the laboratory rest wavelengths for the lines put the 
data into the program linespec to produce a syn- 
thetic template which we call femtemp97 (figure p"4| ). 

Femtemp97 is indeed a better template than emtemp. 
Using the set of 626 duplicate spectra described in 
section 3.1, we can compare the velocity differences 
between pairs directly; the median difference using 
femtemp97 is ~ 9% smaller than in the reduc- 
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Fig. 14. — The emission line template femtemp97. 




Fig. 15. — Residuals in fit to Gaussian error dis- 
tribution for pairs reduced with the femtemp97 tem- 
plate. The contours represent lines of equal residual, 
position in the x,y space represents admixtures of the 
statistical and systematic error models; see text. 



tion using emtemp. The errors for femtemp97 are 
more Gaussian than for emtemp; Figure [l^ shows 
the fit to the two parameter error model (section 
3.2), and may be compared with figure ||. The in- 
terior contour level here is half that in figure Also 
note that the y intercept is about 32km/ s, substan- 
tially less than the 45km/ s in Figure |^ , implying a 
greater reduction in error than the direct compari- 
son of velocity differences would indicate. The me- 
dian ratio (1 + r/ emtemp 9 7 )/(l + r emtemp ) ~ 1.33; 
femtemp97 yields good velocities for a substantial 
number of spectra where emtemp fails. A compari- 



son of 



k f e.mtemp97 
(l - ! - ?" /emtemp97) 



with the # 



8 (1+r) 



error estimator 



shows that 



8 (1+r) 



underestimates the error by 9% 



4.2. absorption line template 

The first absorption line templates used on the 
FAST data were the TD79 vintage ztemp, and the 
NGC4486b template from the MMT spectrograph, as 
well as some secondary MMT templates. Over the 
next year a program of template observations pro- 
vided several high S/N observations of candidate tem- 
plates. These were summed on an object by object 
basis, to form extremely high S/N templates. The 
best of t hes e is the NGC7331 template (fn7331temp) 
(section 3.1). 

We are now in a position to create a better tem- 
plate. Using fn7331temp and femtemp97, we selected 
galaxy spectra where the r value for the reduction 
with fn7331temp was greater than eight, and the r 
value for the reduction with femtemp97 was less than 
3. These are 1959 moderate to high S/N absorption 
line spectra. This set of spectra still contains spec- 
tra of objects with substantial emission. NGC7331 
shows clear emission in the Nil 6583A line, for exam- 
ple. We examined the differences between the veloci- 
ties derived using fn7331temp and femtemp97, where 
(1) the difference was near zero, meaning that there 
was enough emission present to get a correct redshift 
(about 100 spectra) and where (2) the difference was 
near 948/cm/s, meaning that Nil 6583A was confused 
with Ha (about 300 spectra); we removed those spec- 
tra from the sample. 

The remaining 1489 spectra were shifted to a com- 
mon rest velocity, using the fn7331temp velocity; they 
were normalized to the same number of counts; their 
continua were subtracted, using a moderate order 
spline; and fi nally they were summed. The sumspec 
task (section |A.4| ) performed these tasks. The result- 
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ing spectrum had its residual continuum subtracted 
using a high order spline, and was normalized to rep- 
resent the average spectrum. This is the final new 
absorption line FAST template, fabtemp97, shown in 
figure |l6| . 

Fabtemp97 is a better template than fn7331temp. 
The median difference between pairs of spectra from 
the 628 duplicates of section 3.1 is actually ~ 2% 



larger (insignificant) using fabtemp97, because fn7331temp 
is better able to match emission line objects (it has 
Nil visibly in emission, and other lines buried in 
the noise); if one restricts the comparison to pairs 
of spectra where rf ernte mp97 < 3 for each spec- 
trum, fabtemp97 shows a median difference ~ 12% 
smaller than fn7331temp. As with the comparison of 
femtemp97 and emtemp, the differences between the 
new and the old templates are not large. 

Figure |l7] shows the fit to the two parameter er- 
ror model in section ^2 (compare with figure [l(] for 
fn7331temp). The contours in figure [l?] are the same 
as in figure 10; note that even the outer contour does 
not reach the y axis. This suggests that the error 
for these spectra cannot be modeled by a single num- 



ber, k 



5 ^systematic-) 



independent of r; this is not true for 
fn7331temp, or for the emission line templates. A 
comparison of 



error es- 



{l+r fabtemp97 ) 8 (1+r) 

timator shows that | njpp) overestimates the error by 
4%. 

The analyses of the four templates, shown in fig- 
ures g Q H, and all are consistent with a model 
where ~ 20km/ s constant systematic error is com- 
bined with a statistical error determined by the value 
of the r statistic. 
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Fig. 16. — The absorption line template fabtemp97. 




Fig. 17. — Residuals in fit to Gaussian error dis- 
tribution for pairs reduced with the fabtemp97 tem- 
plate. The contours represent lines of equal residual, 
position in the x,y space represents admixtures of the 
statistical and systematic error models; see text. 

We also attempted to make a template for narrow 
lined absorption line objects. We built a template 
by summing, after shifting to a common rest frame, 
several high S/N spectra from a set of M31 globu- 
lar clusters; we call this template fglotemp. Next 
we extracted from the FAST database all spectra 
where rf g i oternp > rf n733 i temp and rf g i oternp > 6 and 
rfemtem P 97 < 3. More than 90% of these spectra were 
calibration stars; we excluded them, and the glob- 
ular clusters themselves, and were left with ~ 300 
galaxy spectra. These we shifted and summed with 
sumspec in the same manner as for fabtemp97, giv- 
ing us a new narrowlined template. This template 
and fabtemp97 were correlated against the 300 spec- 
tra. There was no significant difference in the re- 
sults. From this we conclude that for typical redshift 
survey spectra, from a spectrograph with resolution 
R ~ 1500 there is no need to have a narrow lined 
template. 

4.3. velocity zero point 

We have developed a new methodology for defin- 
ing the velocity zero point. Previous methods have 
depended on external calibrators, such as 21 cm mea- 
surements; our new methods are fully internal, and 
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Fig. 18. — Spectral type differences in determining 
velocity zero points. The top panel shows velocities 
for 75 observations of NGC4486b, the solid histogram 
using the fn7331temp template, and the dotted his- 
togram using the ztemp template. The x axis has 
been shifted so that the median redshift obtained us- 
ing ztemp is zero. The bottom panel shows the same 
information for 116 different spectra of M31. No shift- 
ing of zero points can make both sets of histograms 
agree. 



should minimize velocity offsets due to spectral type 
differences. 

Figure |l8| illustrates the problem with spectral type 



diff erence. With the FAST spectrograph we observed et al., 199Cj ) 



M31 on 116 occasions, and NGC 4486b on 75 occa- 
sions. Each of these spectra was reduced using ztemp 
and fn7331temp. The resulting velocities were shifted 
so that, for each galaxy, the median redshift obtained 
by using the ztemp template was zero. The bottom 
panel of figure [l^ shows the results for M31; the solid 
histogram shows the result of correlating the 116 spec- 
tra with fn7331temp, the dotted histogram shows the 
result of correlating these spectra with ztemp. Simi- 
larly, the top panel of figure 18 shows the results for 
NGC 4486b, with the solid histogram representing 75 
fn7331temp correlations, and the dotted histogram 75 
ztemp correlations. 

No change in the zero point of either template can 



make both sets of histograms agree; any setting of 
the zero point by matching velocities for one object 
would make the systematic difference between differ- 
ent template reductions for the other object worse. 

We choose to define the zero point of our velocity 
system to minimize the systematic difference between 
the two main types of galaxy spectra: emission line 
spectra and absorption line spectra. To insure that 
our internal procedure matches the "true" system we 
need only make the reasonable assumption that we ac- 
curately know the rest velocities of the main spectral 
lines in galaxies, such as Ha and [O III]. We must 
also be certain that the internal wavelength system 
of the spectrograph matches the external wavelength 
system of the sky (section ||). 

Our basic procedure is to force the median differ- 
ence between the emission line velocity and the ab- 
sorption line velocity, for those galaxies which strongly 
show both sets of features, to zero. Because the ab- 
sorption line velocity comes from the K giant stars in 
the bulge, primarily, while the emission line velocity 
comes mainly from HII regions in the disk, which are 
moving at ~ ±200km/s with respect to the central 
bulge, there is no reason to expect that the absorption 
line velocity and the emission line velocity should be 
the same for any particular object. These differences 
should be randomly distributed about zero, however. 

We began by setting the velocity for the fn7331temp 
template so that the median velocity difference be- 
tween the fn7331temp velocity and the femtemp97 ve- 
locity, for spectra where r/ emtemp9 7 > 5 andr/ n 733it emp > 
5 was identically zero. This yields a velocity of 
797km/ s for NGC 7331, which may be compared with 
820±3/cm/s from the 21 cm observations ( BottineflT] 



As described in section |4^, fabtemp97 was created 
by shifting 1489 spectra to the rest frame defined by 
their individual correlations with fn7331temp, then 
summing them. A comparison of fabtemp97 velocities 
with emission line velocities should give an indication 
of the stability of the zero point technique. 



Figure 19 shows the difference in velocities between 
fabtemp97 and femtemp97 as a function of observa- 
tion date, for 1787 FAST spectra where rf em tem P 97 > 
5 and rf a btem P 97 > 5. The median difference is 
7.1km/ s, with an interquartile range of 62km/ s. Se- 
lecting subsets of the spectra with higher S/N ratios 
has no significant effect on this result. This difference 
is ~ 0.1 pixel, and is probably due to systematic, 
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Fig. 19. — The difference between velocities obtained 
with the femtemp97 template and the fabtemp97 tem- 
plate, for 1787 FAST spectra where both reductions 
had r > 5, as a function of observation date. 



Fig. 21. — The difference between velocities ob- 
tained by fitting O[III]5007A and correlating with 
the femtemp97 template, for 3527 spectra where r > 5 
and 0[III] was fit, as a function of observation date. 
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Fig. 20. — The difference between velocities obtained 
by fitting Ha and correlating with the fabtemp97 tem- 
plate, for 1514 spectra where r > 5 and Ha was fit, 
as a function of observation date. 
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Fig. 22. — The difference between velocities obtained 
by fitting Ha and correlating with the femtemp97 
template, for 4833 spectra where r > 5 and Ha was 
fit, as a function of observation date. 
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Fig. 23. — The difference between velocities obtained 
by fitting 0[III] 50071 and fitting Ha, for 2008 spec- 
tra where Ha and 0[III] were both fit, as a function 
of observation date. 



non-linear errors in our wavelength scale. Figure |2Q 
shows the difference between the fabtemp97 velocity 
and the velocity found by fitting Ha in emsao for a 
representative subset of the data. The median dif- 
ference is — 2.1 km I s, consistent with zero. Figure p2| 
shows the difference between the femtemp97 veloc- 
ity and the Ha velocity, — 9. 51 /era/ s, and figure ^l| 
shows the difference between the femtemp97 velocity 
and the [OIII] 5007Avelocity, -0.16km/ s. These dif- 
ferences are very similar to the median differences be- 
tween individual lines measured in the same spectrum 
with emsao, Figure [23] shows the difference between 
fits to Ha and fits to 0[III], the median difference is 
QAkm/s. 

We conclude that our construction of fabtemp97 
gives a velocity zero point equal to the velocity zero 
point of the emission line system; this equivalence is 
as accurate as our ability to define the wavelength 
system using standard lamps and polynomial fits. 

5. Testing spectrograph zero point and sta- 
bility using emsao and night sky spectra 

Before 1995, we measured redshifts for emission 
line spectra using emsao. Now we cross-correlate 

Odnst an emission line template with xcsao (section 
). emsao is used to check for errors in emission- 
line correlations, to interactively reduce emission- line 
spectra, to automatically measure equivalent widths, 
line heights, and line widths for sets of emission line 



spectra, and to do various custom projects using al- 
ternate line lists (e.g. to study QSOs). 

We used emsao to calibrate the Z-Machine and 
FAST spectrographs for stability and zero point by 
measuring the apparent velocity of the night sky 
lines. These are from forbidden oxygen air glow lines 
from the upper atmosphere and mercury and sodium 
emission from artificial sources such as street lights. 
Sodium is a blended doublet and we do not know 
the effective wavelength a priori] so we cannot use it 
to establish a zero point. Using the zero point de- 
fined by 0[I] 5577A we measure an effective wave- 
length of 5891. 2 A using FAST, with the Z-Machine 
data in agreement. We also obtain a systematic dif- 
ference between the oxygen airglow lines and the mer- 
cury streetlamp lines of about 20km /s with both 
spectrographs. We assume that this effect is due to 
differences between the effective wavelengths of Hg 
in street lamps and calibration lamps. We therefore 
choose to calibrate our zero points with the oxygen 
airglow lines. 

Figure 24 shows the results for the Z-Machine. The 
top panel is a typical sky spectrum, the middle panel 
shows the apparent velocity of the 0[I] 5577A line, 
and the bottom panel shows the apparent velocity of 
0[I] 6300 A . The velocities are shown as a function 
of observation date, over 15 years. Each point is the 
median of a single nights observations; nights with 
fewer than 10 observations were excluded. 

The night sky emission lines were routinely used 
to supplement the HeNeAr calibration lamp lines in 
the Z-Machine reductions, and a seventh order poly- 
nomial was used to fit the lines. The 0[I] 5577A line 
does not show the same general pattern as the 0[I] 
6300A line. All the other lines (Hg 5461, NaD, 0[I] 
6363) do show the same behavior as the 0[I] 6300 
line. 0[I] 5577A is in a portion of the spectrum with 
no strong HeNeAr calibration lines, which, along with 
the high order polynomial, locally matches the cali- 
bration to the position of 0[I], instead of the system 
defined by the calibration lamp. 

While the 0[I] line at 6300 A was also used in the 
wavelength calibration, it is in a region with many 
strong Ne lines, and should have negligible effect on 
that calibration. The apparent velocity of this line 
should be a good measure of the difference between 
the instrumental zero point and the true sky. The 
scatter and long-term changes in this line measure 
the stability of the instrument and data reduction 
procedures. Many of the large jumps in apparent ve- 
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Fig. 24. — emsao reduction of night sky lines from 
the Z-Machine. The top panel shows a typical night 
sky spectrum, with the lines marked. The middle 
panel shows the apparent velocity of 0[I] 5577A as a 
function of observation date. Each point is the nightly 
median, where nights with fewer than 10 observations 
are excluded. The bottom panel is similar to the mid- 
dle, but shows the apparent velocity of the 0[I] 6300A 
line. 
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Fig. 25. — emsao reduction of night sky lines from 
FAST. The top panel shows a typical night sky spec- 
trum, with the lines marked. The middle panel shows 
the apparent velocity of 0[I] 5577A as a function of 
observation date. Each point is the nightly median, 
where nights with fewer than 10 observations are ex- 
cluded. The bottom panel is similar to the middle, 
but shows the apparent velocity of the 0[I] 6300 A 
line. Note that the scale of the ordinate is half that 
in figure \24 
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locity can be attributed to known hardware changes; 
the largest feature in figure |3, the 5.5 year slow in- 
crease in velocity of 0[I] 6300A from ~ — 30km/ s to 
~ +10/cra/s is of unknown origin. 

Figure ^5| shows the results for FAST. Again the 
top is a typical sky spectrum, the middle shows the 
apparent velocity of the 0[I] 5577A line, and the bot- 
tom the apparent velocity of the 0[I] 6300A line. The 
abscissa covers four years, and the scale of the ordi- 
nate is half that of figure |24|. 

The night sky lines are not used in the wavelength 
calibration of FAST, and the HeNeArFe calibration 
lamp lines are fit with a third order polynomial. Both 
0[I] 5577A and 0[I] 6300A should be good measures 
of the instrumental zero point and the stability of the 
instrument and reduction techniques. 

Both lines show essentially the same behavior, es- 
pecially since the CCD was changed in September 
1994. The scatter in these line positions is less than 
20% of the scatter for the same lines in the Z-Machine. 
Over the four years of FAST operation an additive off- 
set of 7.5±2/cm/s brings the instrumental system into 
agreement with the true system of the sky. 

FAST is a remarkably stable instrument. The scat- 
ter in theO[I] 5577A line can provide a good measure 
of the errors due to a combination of instrumental in- 
stability, wavelength calibration, and line fitting. We 
measure the error in fitting a single line by dividing 
the scatter in the apparent velocity difference between 
0[I] 6300 A and 0[I] 6363 A by \/2. Subtracting it in 
quadrature from the scatter in 0[I] 5577A , yields the 
error due to the interaction of instrumental instabil- 
ity with the wavelength calibration. For the entire 
four year period this is 2km/ s, where most of the 
error comes from systematic changes associated with 
changing the CCD, changing the linelist for the He- 
NeArFe calibration, and changing the dewer. For the 
most recent 18 months, since the dewer change, the 
error from instrumental unstability is unmeasurably 
small, is consistent with zero, and has a 2a upper 
limit, by bootstrap resampling, of 0.7km/ s. 

6. Other Methods and Comments 

Press (1995) has suggested a new methodology 
for determining redshifts. One first reduces a set of 
galaxy spectra (shifted to a common rest frame) to a 
set of orthogonal basis vectors, using singular value 
decomposition, SVD. Next, using the fast numerical 
methods of Rybicki and Press (1995) the most signif- 



icant of these vectors are repeatedly fit to a spectrum 
with unknown redshift, with each fit being at a differ- 
ent redshift. For each fit x 2 is calculated; the redshift 
corresponding to a minimum % 2 is the redshift of the 
galaxy. 

Recently Glazebrook, et al (1998) have developed a 
nearly identical scheme; the difference is that rather 
than calculate \ 2 exactly, they use the simplifying 
assumption that the correlation function may sub- 
stitute for x 2 • Thus they can treat the (SVD de- 
rived) eigenvectors as templates in a cross correla- 
tion program (such as xcsao), and obtain a final red- 
shift by summing the correlation functions in quadra- 
ture, weighted by the eigenvalues of each eigenvector- 
template. 

Both Press (1995) and Glazebrook, et al (1998) 
have suggested that the coefficients of the fits can be 
used to classify the spectra, and that the position of a 
spectrum in coefficient space may be used to develop 
confidence measures in the derived redshift. Recently 
Bromley, et al (1997) used this technique to classify 
spectra from the LCRS (Shectman, et al, 1996). 

We have not adopted the SVD method for creating 
templates, nor the x 2 minimization technique for de- 
termining redshifts. We have, following Press's (1995) 
suggestion, developed methods to use best fit param- 
eters for rough classification and blunder discovery. 

We do not implement the Glazebrook, et al (1998) 
method for a number of reasons. Their approximation 
assumes the variance in a spectrum is not a function of 
wavelength. This assumption is clearly incorrect, and 
is one of the reasons why continuum subtraction is 
superior to continuum division in the low S/N regime 
(section ^ ) . The use of the eigenvalue as a weight in 
the combination of correlation functions overweights 
emission line correlations when reducing absorption 
line spectra, thus making them even more susceptible 
to shot noise. But finally, as the exact methods are 
available (Press, 1995), we see no advantage in the 
approximation. 

Press's (1995) x 2 minimization technique, enabled 
by the Rybicki and Press (1995) algorithm, shows 
substantial promise. Especially for the case where 
the rest frame of the template is much different from 
that of the unknown galaxy we expect that a rigor- 
ous accounting of observational errors as a function 
of wavelength will be important; the new, deep red- 
shift surveys of the next decade, e.g. with Hectospec 
(Fabricant, et al, 1994), will test whether the x 2 tech- 
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niques will obtain better results for high redshift ob- 
jects. We do not now implement the x 2 minimization 
because we believe that current redshift survey data, 
with z ~ 0.2 would not benefit; some confirmation of 
this view comes from the preliminary results of Press 
(1995, 1997), who finds no improved ability to obtain 
redshifts from low S/N spectra in the LCRS (Sheet- 
man, et al 1996) compared with the original reduc- 
tion, done by H. Lin using RVSAO 1.0. 

We believe our template creation techniques are su- 
perior to a SVD decomposition of a group of spectra 
for underlying physical reasons. Emission line spectra 
and absorption line spectra arise from independent 
physical causes in different physical locations within 
a galaxy; basis vectors which are admixtures of ab- 
sorption and emission line spectra make no physical 
sense. Additionally emission line velocities and ab- 
sorption line velocities are not identical for any par- 
ticular galaxy, as a perusal of optical rotation curves 
(e.g. Barton, et al, 1998) clearly shows. 

The simultaneous use of more than two templates 
or basis vectors to obtain redshifts is unnecessary. We 
performed a SVD decomposition on the 1489 pure 
absorption line spectra used to create fabtemp97; the 
first eigenvector was essentially identical to fabtemp97, 
the next three represented small differences in the 
continuum subtractions, the fifth eigenvector shows 
the H+K lines, which are systematically weaker in 
the early data, before a blue sensitive chip was in- 
stalled, the next couple of eigenvectors also represent 
continuum differences. The systematic residuals from 
the reduction dominate the higher dimensions of the 
SVD decomposition. 

As a practical matter the limiting factor in obtain- 
ing redshifts for faint objects is the ability to get red- 
shifts for weak lined absorption line systems. Essen- 
tially, if one observes long enough with a fiber instru- 
ment to get redshifts for the absorption line spectra, 
the spectra with emission lines will all yield redshifts 
with almost any technique. 

We therefore will use fabtemp97 and femtemp97 as 
our fitting functions, we fit each spectrum as a linear 
combination of these two templates. We take as a 
representative sample the last 2000 spectra observed 
as part of the 15R survey (Geller, et al, 1998). Fig- 
ure ^6| shows the absorption line vs. emission line fit 
coefficients. The solid dots are for r/ em < 3 and the 
open circles are for r/ em > 3. There is a clear locus, 
outliers could be easily discovered and removed from 
any fully automatic data reduction process. Com- 
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Fig. 26. — The coefficient of the fit to fabtemp97 
vs the coefficient of the fit to femtemp97, for 2000 
spectra from the 15R survey (Geller, et al (1998)). 
About a dozen objects have coefficients not in this 
range, and would be looked at manually. The solid 
dots are spectra where rf emtemv < 3, the open circles 
are for spectra where r/ em£emp > 3. The diagram 
can be viewed as a crude absorption line strength vs 
emission line strength classification diagram. 

bined with more traditional measures, e.g. equivalent 
widths and line ratios, figure |2^ can form the basis 
for a spectroscopic classification. The x and y axes 
roughly measure absorption line strength (metallicity 
perhaps) and emission line strength. 

Figure |2^ shows the relation between the fit resid- 
uals and the inverse of the 1 + Tj a h statistic for the 
fabtemp97 reduction. The symbols have the same 
meaning as in figure [26| There is a very good corre- 
lation of fit residuals with 1/(1 + r) for spectra where 
ffem < 3; for objects with emission lines the correla- 
tion is much weaker. 

Figure |2^ shows the relation between the fit resid- 
uals and the inverse of the 1 + r/ em statistic for the 
femtemp97 reduction, the symbols here are dots if 
ffab > 3 and open circles if rf a b < 3. The correla- 
tion of fit residuals with 1/(1 + r) for spectra where 
ffab < 3 is weak, and for objects with absorption line 
velocities nonexistent. 

These two plots, along with the figures Il5| and Il7| in 
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Fig. 27. — The relation of the residuals to the fit to 
fabtemp97 to the inverse of l+rf a btemp97 The symbols 
have the same meaning as in figure pq. 
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Fig. 28. — The relation of the residuals to the fit 
to femtemp97 to the inverse of 1 + rf ernte m P 97- The 
solid dots are for spectra where rf a btem P 97 > 3 and 
the open circles are for spectra where rf a btem P 97 < 3. 



section which show the statistical and systematic 
error models, demonstrate that for typical redshift 
survey spectra, residuals for absorption line objects 
are mainly determined by signal to noise; for emission 
line objects, the residuals are mainly systematic. 

7. Conclusion 

We have demonstrated the techniques which we 
have developed in the RVSAO suite to create a sys- 
tem for the accurate, automated reduction of spec- 
tra for galaxy redshifts and stellar radial velocities. 
More than half of all published redshifts have been 
measured using these techniques, as well as a large 
number of stellar radial velocities. 

The correlation method for obtaining redshifts can 
be successfully extended from absorption line spectra 
to emission line spectra, with a substantial improve- 
ment in effectiveness over the previous method for ob- 
taining emission line redshifts, automated line fitting. 
The reduction of emission-line spectra requires differ- 
ent reduction steps than absorption line correlations. 
Emission line correlation redshifts are susceptible to 
blunders due to the presence of cosmic rays. However, 
using automated line fitting (emsao) and absorption 
line correlation velocities the blunder rate can be kept 
near zero, with the degree of automation kept high. 

We have developed new techniques for calibrating 
and characterizing the blunder rate and the individ- 
ual errors in redshift measurements. The blunder rate 
for RVSAO reductions can be kept near zero by the 
use of some simple heuristics to identify possible mis- 
takes. For typical redshift survey data from the FAST 
spectrograph the automation rate is 95%. Our self 
calibrating internal error estimator is accurate to ~ 
20%. Large, stable surveys enable development of 
more accurate and stable error estimators. 

We have developed new methods for creating, cal- 
ibrating, and using galaxy redshift templates. We 
have created an emission line template, femtemp97, 
having the median properties of a large set of strong 
emission line spectra. We have created an absorption 
line template, fabtemp97, having the mean proper- 
ties of a large set of absorption line spectra showing 
no sign of emission. These spectra arise from phys- 
ically distinct processes, and can be used to form a 
pair of basis vectors to perform a 2-D spectral clas- 
sification. We have developed a new method for es- 
tablishing the zero point for redshift observations, a 
method which minimizes the systematic differences 
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between emission and absorption line redshifts. This 
zero point is determined as accurately as we can es- 
tablish the wavelength calibration using standard He- 
NeArFe lamps. We have shown a technique for mea- 
suring and eliminating differences between the instru- 
mental zero point and the true zero point. 

We have shown improved techniques for a number 
of the substeps necessary to obtain accurate redshifts, 
including: removal of emission(absorption) lines when 
correlating against an absorption(emission) line tem- 
plate; suppression of the night sky lines; supression of 
the continuum; design of the Fourier filter; and zero 
padding of spectra. 

The rapid development of large aperture telescopes 
with multi-object spectrographs presents substantial 
challenges for redshift and radial velocity reductions. 
Reducing one or two orders of magnitude more spec- 
tra of objects one or two orders of magnitude fainter 
while maintaining high quality control standards and 
minimal personnel costs is clearly a difficult problem. 
RVSAO provides a solid methodological and software 
basis to meet these new challenges. 
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A. The Elements of RVSAO 

The RVSAO package consists of six IRAF tasks 
xcsao, emsao, linespec, sumspec, contpars, and 
bcvcorr. Each of these tasks is controlled by a set of 
user settable parameters. In this appendix we list and 
describe all the relevant parameters for these tasks, 
and demonstrate their use in a series of examples. 



A.l. Cross-Correlating a Spectrum in xcsao 

Digital cross-correlations in the RVSAO system are 
performed by xcsao, basically following the prescrip- 
tion of TD79, with a large number of refinements and 
additions, xcsao is capable of reducing a wide variety 
of input spectra, but must have its many parameters 
properly set. Param eters for xcsao are set in its pa- 
rameter list (Figures |A29| and |A30), in the contpars 
program's parameter list (Figure A67 ), and by spe- 
cial instructions to the program in the headers of the 
template spectra (Figure |A3l| ). In this section we 
examine in detail the elements of xcsao, and we de- 
scribe how the parameters are used in the code, along 
with the algorithmic details of their use. We follow 
the cross-correlation of a single spectrum against two 
very different templates, one showing primarily ab- 
sorption features and one showing only emission fea- 
tures, to demonstrate the ways in which processing 
varies in response to the template spectra being used. 

For each object spectrum file from the input list 
spectra, and/or each aperture specified in the aperture 
list specnum, the fitting subroutine is called. Spec- 
trum files are all read from the directory specified by 
specdir, but full or relative pathnames may be used in 
spectra, and specdir may be null. Spectra which are 
in flux units instead of counts should be renormalized 
by setting renormalize to yes. If obj-pl ot is yes, the 
object spectrum is plotted as in Figure |A32| , and the 
plot is kept on the screen available for zooming and 
editing until a "q" is typed. If fixbad is "yes," regions 
specified in the file named by badlines are replaced by 
straight lines connecting the adjacent pixels, and the 
image is plotted again if objjplot is "yes." 

Band tempband of each template spectrum in the 
list templates and list of multispec apertures, temp- 
num, is loaded. Template files are all read from the di- 
rectory specified by tempdir, but full or relative path- 
names may be used in templates, and tempdir may be 
null. If echelle is yes, tempnum is ignored, and the 
multispec lines used for templates track those used 
for object spectra. If tempjplot is yes, the template 
spectrum is plotted. 

A zero-point redshift is computed by adding the 
solar system barycentric velocity correction, from a 
source specified by sveLcorr, the redshift of the tem- 
plate from the VELOCITY parameter in the tem- 
plate header, an optional template-dependent veloc- 
ity shift from the TSHIFT parameter in the header, 
and an optional constant velocity shift from the tshift 
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fabtemp97, femtemp97 



archive$templates / 



badlines . dat 



(s_abs_r= 



4096 
tempfile 



INDEF 
INDEF 
ablines . dat 
emlines . dat 
rvsao$lib 



List of input spectra 

i aperture range if multispec or 2-D file 
Spectrum band if multispec file 
Directory for input spectra 
Cross-correlate spectrum (yes or no) 
List of template spectra 

Template aperture range if 2-D or multispec file 
Template band if template is multispec file 
Directory for template spectra 

Echelle spectra (template tracks spectrum) (yes or no) 

Starting wave length to correlate 

Ending wave length to correlate 

Plot the raw object data 

Plot the filtered correlation function 

Write the filtered correlation function to a file 

Eliminate portions of spectrum with bad lines (yes or no) 

Lines to eliminate from spectrum 

Remove emission or absorption lines (yes, no, tempfile) 
Remove emission or absorption lines (yes, no, tempfile) 
Spectrum absorption line rejection in sigma of fit 
Spectrum emission line rejection in sigma of fit 
Template absorption line rejection in sigma of fit 
Template emission line rejection in sigma of fit 
Fraction of spectrum to taper 

Renormalize spectrum before transform (yes or no) 
Number of cols into which to rebin 
Rebin interpolation mode 

Pad transforms with zeroes to lower noise (yes, no, tempfile) 
Bin number of low frequency cutoff 

Bin number where high frequency starts to slope off 
Bin number of high frequency cutoff 

Initial velocity source (zero guess zguess corr em comb) 
Initial velocity guess in km/sec or z 

Number of times template is shifted to match object velocity 
night to night zero point velocity shift 

Spectrum velocity correction (none, file, helio, bary, hf ile) 
Template velocity correction (none, file, helio, bary, hfile) 
Peak-fitting mode (1-parabola 2-quartic 3-cosx/l+x A 2) 
Fraction of peak or number of points for peak fitting (0=FWHM) 
Half-width of search region for graphically set peak 
Minimum allowable velocity in km/sec 
Maximum allowable velocity in km/sec 
Report mode (l=normal 2=one-line) 
List of log files 

Save result in I RAF data file header (yes or no) 
Save results in binary archive record (yes or no) 
Smooth displayed spectrum (0=no else number of times) 
Center velocity of output graph 
Velocity half-width of output graph 
Absorption line list 
Emission line list 
Directory for line lists 

Template absorption line rejection in sigma of fit 
Template emission line rejection in sigma of fit 



Fig. A29. — Parameter list for xcsao. 



Keyword Value 
CHOPEM T or F 



Usage 

If T and xcsao.s_emchop=tempfile, replace all pixels in the object spectrum 
which are less than contpars.s_abs_reject or greater than contpars.s_em_reject 
standard deviations from the continuum with the value of the continuum fit at that 
wavelength. If T and xcsao.t_emchop=tempfile, replace all pixels in the template 
spectrum which are less than contpars.t_abs_reject sigma or greater than 
contpars.t_em_reject standard deviations from the continuum with the value of 
the continuum fit at that wavelength. 



#sig If CHOPEM is true and xcsao. s_emchop=tempfile, this limit is _ 

minimum for the object spectrum instead of contpars.s_abs_rejec 



used as the 



If CHOPEM is true and xcsao. s_emchop=tempfile, this limit is used as the 
maximum instead of contpars.s_em_reject. It does not override 
contpars.t_em_reject for the template spectrum. 



SUBCONT TorF 



DIVCONT T or F 



ZEROPAD TorF 



This flag determines how the Fourier transforms are filtered: 

or not present = Filter both the template and the object spectrum transforms 

using xcsao.low_bin, xcsao. top_low, xcsao.top_nrun, and xcsao.nrun to define a 
cosine-bell. 

1 = Do not filter the template spectrum transform, but filter the object spectrum transform 

2 = Do not filter the high frequencies of the object and template spectrum transforms. 

3 = Do not filter the high frequencies of the object spectrum transform and do not 

filter the template spectrum transform at all. 

Override the pkfrac parameter for just this template by fitting the correlation peak 
from its maximum value down to this fraction of the maximum. The most common 
use of this parameter is to fit only the top quarter (set to ) of a correlation with an 
emission line template to avoid fitting side lobes caused by close emission lines 
in the template, such as Nl and N2 next to H-alpha. 

If true or not present, fit and remove the continuum from the template; if false, do 
not fit or remove continuum from template. This flag does not affect removal of 
the continuum from the object spectrum. 

If T, divide, instead of subtract, the continuum from template and spectrum. 

If the ZEROPAD parameter in XCSAO is set to tempfile instead of yes or no, and 
this template spectrum header parameter is T, the spectrum and template 
transforms are both padded with an equal number of zeroes before being 
cross-correlated. If it is F, the transforms are not zero-padded. 



Fig. A31. — These template spectrum header pa- 
rameters can be used to control the cross-correlation 
process to correlate object spectra against both ab- 
sorption line templates and emission line templates 



m one xcsao run. 



4096 
tempfile 



tempfile 



Fraction of spectrum to taper 
Renormalize spectrum before transform {yes 
Number of cols into which to rebin 
Rebin interpolation mode 

Pad transforms with zeroes to lower noise { 
Bin number of low frequency cutoff 
Bin number where low frequency reaches one 
Bin number where high frequency starts to slope off 
Bin number of high frequency cutoff 

Initial velocity source {zero guess zguess corr em comb) 
Initial velocity guess in km/sec or z 

Number of times template is shifted to match object velocity 
Spectrum velocity correction (none, file, helio, bary, hfile) 
Template velocity correction (none, file, helio, bary, hfile) 
Peak-fitting mode (1-parabola 2-quartic 3-cosx/l+x A 2) 
Fraction of peak or number of points for peak fitting (0=FWHM) 
Half-width of search region for graphically set peak 
Minimum allowable velocity in km/sec 
Maximum allowable velocity in km/sec 
Report mode (l=normal 2=one-line) 
List of log files 

Save result in IRAF data file header (yes or no) 

Save results in binary archive record (yes or no) 

Center velocity of output graph 

Velocity half-width of output graph 

Absorption line list 

Emission line list 

Directory for line lists 

Display modes (l=spectrum+xcor 2=spectrum+lines) 
Display results on <device> (yes or no) 
Display device 

Wait for keyboard response after each object (yes or no) 

Automatic hardcopy of results on <plotter> (yes or no) 

Hardcopy output device 

Plot the template spectra 

Plot the continuum-subtracted data 

Plot the apodized data 

Plot the fft of the data 

Plot the fft of the fft to show what is be: 
Plot the unfiltered cross-correlation 
Display intermediate results 
Graphics cursor input 



Fig. A30. — Parameter list for xcsao continued. 




Wavelength in Angstroms 



Fig. A32. — Object spectrum showing both emission 
and absorption lines as plotted by xcsao if obj.plot is 
yes. 
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parameter. The template spectrum's barycentric ve- 
locity correction, from a source specified by tveLcorr, 
is subtracted because the template spectrum's ob- 
served velocity, not its corrected one, gives that spec- 
trum the redshift which we are comparing. An ini- 
tial redshift source may be specified by veLinit; if 
"guess" , this is from czguess; if "zero" the initial ve- 
locity is 0; otherwise it can be read from an object 
spectrum header parameter VELOCITY ("combina- 
tion"), CZXC ("correlation"), or CZEM ("emission"). 
If such an initial redshift has been called for, or if this 
is the second pass or greater (nzpass > 1), the tem- 
plate log- wavelength limits are shifted by that initial 
redshift (on the first pass) or the current correlation 
redshift. The wavelength region over which the tem- 
plate and object spectra overlap is computed. If the 
wavelength in Angstroms specified by stJambda, is 
greater than the blue limit of the overlap region, it be- 
comes the new limit. If the wavelength in Angstroms 
specified by endJambda, is less than the red limit of 
the overlap region, it becomes the new limit. The 
two spectra are rebinned into log-wavelength with a 
number of pixels, set by ncols using an interpolation 
mode specified by interpjmode. 

First, the continuum, and, optionally, emission 
and/or emission lines, are removed from the rebinned 
object and template spectra. Parameters for fitting 
the continuum are in the IRAF pset task named con- 
tpars (see section |A.5| ). Emission and/or absorption 
lines may be removed from either or both of each 
object-template pair, s.emchop controls whether lines 
will be removed from the object spectrum. If the 
template spectrum header parameter SUBCONT is 
present, its value overrides that of S-emchop. s.absrej 
and S-emrej set the lower and upper acceptable limits 
for object spectrum pixels to be used in the continuum 
fit. If lines are rejected, as they are when the object 
spectrum we are following is to be correlated against 
an absorption line temp late, the rejected data points 
are plotted as in figure A33 if contsub-plot is yes. A 
graph of th e continuum-removed data, as shown in 
Figure |A34 is displayed if the contsub-plot parame- 
ter is yes. 

Template pixels with values outside of the lower 
t-absrej and upper t-emrej acceptable limits, in stan- 
dard deviations from the continuum fit, are replaced 
by continuum values in the template spectrum if 
t-emchop is set to yes. The continuum is then sub- 
tracted from the template spectra. 

Template and object spectra are then apodized, 




Wavelength in Angstroms 



Fig. A33. — Emission lines removed from the object 
spectrum before cross-correlation against an absorp- 
tion line template spectrum in xcsao. 




5000 5500 
l in Angstroms 



b) 




Fig. A34. — The object spectrum with a) its contin- 
uum subtracted and emission lines removed for cross- 
correlation against an absorption line template and b) 
its continuum subtracted, but its emission lines kept 
for cross-correlation against an emission line template 
spectrum in xcsao. 
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a) 



b) 




Fig. A35. — The object spectrum apodized at 
each end, ready to be Fourier transformed for cross- 
correlation against a) absorption and b) emission tem- 
plates in xcsao. 



a) 




b) 




Fig. A36. — The object spectrum Fourier- 
transformed over the wavelength ranges of the a) ab- 
sorption and b) emission template spectra in xcsao. 



tapered at each end for the fraction belLwindow of 
the entire spectrum. If zeropad is yes, both spectra are 
padded with an equal length of null (zero) spectrum. 
A graph of the continuum-removed, apod ized object 
spectrum is displayed, as shown in Figure |A35| , if the 
apodize-plot parameter is yes. 

The spectra are then Fourier transformed. The 
Fourier power spectra of the obje ct an d te mpla te 
spectra are displayed as in Figures A36 and A37 if 
ffljplot is set to yes. 

The transformed spectra are then filtered with a 
cosine-bell filter. The low frequencies are filtered from 
low-bin to topJow and the high frequencies are filtered 
from topjnrun to nrun. The template header param- 
eter FI-FLAG controls whether the template trans- 
form is filtered and whether the high-frequency filter 
is turned off for both template and object transforms 
to leave in emission lines. If tfftjplot is set to yes, 
the filtered Fourier power spectra of the o bject and 
template spectra are displayed as in Figures |A38| and 

The filtered transforms are then cross-correlated 



and normalized. If uxcorjplot is yes, the unfiltered 
correlation is displayed. If xcorjplot is yes, this re- 
sult is displayed as in Figure |A4Q| , and a specific peak 
may be selected using the cursor if curmode is yes. 
In that case, the maximum value within pksrch pixels 
of the cursor position is used. Otherwise, the high- 
est correlation peak between the velocities minvel and 
maxvel is used. The redshift is calculated by fitting 
a parabola or similar function specified by pkmode to 
the portion of the peak above pkfrac of the maximum 
value of that peak. The R-value and error are com- 
puted, and control returns to XCFIT to set up the 
template for the next pass. 

After all of the template spectra have been corre- 
lated against an object spectrum, the template with 
the highest R-value is selected. The results are dis- 
played as text to the devices specified by lo gfiles in the 
format specified by report-mode. Figure A41 shows 
the default report. 

If displot is yes, the object spectrum and, option- 
ally, the selected correlation peak, are plotted to de- 
vice in the format specified by dispmode. Figure A42 
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b) 




b) 




Fig. A37. — Fourier transforms of a) absorption and 
b) emission template spectra, in xcsao. 



shows the summary graph if dispmode is 1. Figure 
|A43| shows the dispmode 2 summary graph, with emis- 
sion and absorption lines labeled. If hardcopy is yes, 
the same graph is sent a printer. 

If nsmooth > 0, the object spectrum is smoothed 
by a 1-2-1 sliding filter nsmooth times for display pur- 
poses only. This smoothing may be changed interac- 
tively using the g command in cursor mode. The 
filtered cross-correlation with the best R- value is dis- 
played centered on the redshift cvel (in km/sec) with 
a width in km/sec of dvel. If cvel is INDEF, the fit 
redshift is used; if dvel is INDEF, the width is set to 
20 times the peak width. If the correlation is not dis- 
played, absorption lines {ablines— yes) and/or emis- 
sion lines (emlines— yes) may be labeled fro m line lists 
in the directory linedir, as shown in Figure |A43 . 



If curmode is yes, the user can interact with the dis- 
play using the terminal cursor to zoom in on portions 
of the spectrum, rerun the cross-correlation, change 
the display format, edit the spe ctrum , or several other 
functions. For example, figure A44 shows the corre- 
lation result for the second best template, selected 
using the T cursor command. 



Fig. A38. — The object spectrum transform filtered 
for correlation against a) absorption and b) emission 
templates in xcsao. Removing high frequencies ad- 
versely affects the shape of narrow emission lines, so 
no high frequency filtering is done for emission line 
correlations 

If savejuel is yes, cross-correlation redshift results 
are written into the object spectrum image header 
in a form appropriate to the spectrum format: two 
entries plus one per template if multispec; otherwise 
one value per keyword, as in Figure |A45| . 

A. 2. Fitting Redshifted Lines in a Spectrum 
in emsao 

While xcsao can now use emission line templates, 
it is still useful to measure redshifts of emission line 
spectra directly. For large surveys, the interactive de- 
termination of line centers and calculation of redshifts 
by a program like IRAF's splot, is simply too slow, 
emsao, a companion to the cross-correlation task xc- 
sao, was written to find emission lines automatically, 
compute redshifts for each identified line, and com- 
bine them into a single radial velocity. The results 
may be graphically displayed or printed, saved to a 
file, and/or stored in the spectrum file header, emsao 
is designed to run with minimal human intervention, 
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b) 




Fig. A39. — Filtered Fourier transforms of a) absorp- 
tion and b) emission templates in xcsao. No high 
frequencies are filtered out of the emission template 
transform. 



b) 




Fig. A40. — Normalized, filtered cross correlations of 
the Fourier transforms of the object spectrum against 
the Fourier transforms of a) absorption and b) emis- 
sion templates in xcsao. 



1995. 1122. 0113. 311. 062165. ms Object: 311.062165 

RA: 08:13:29.34 Dec: 29:02:12.0 $ . If 23-Nov-1995 10:37:37.00 =HJD2450044 . 9457 BCV 
: 25.30 

3679. 9A- 7551. 7A 4096x2 points, filter: 5 20 125 250 apodize 0.05 fit 0.5 best 2 
Temp: fabtemp97 vel : 0.00 tsh: 0.00 HCV: 0.00 Peak: 14415.2 h:0.22 

3 R: 2.88 CZ: 14441.7 +/- 71.592 1 +tc 

Temp: femtemp97 vel: 0.00 tsh: 0.00 HCV: 0.00 Peak: 14369.4 h:0.89 

1 R: 16.38 CZ: 14395.9 +/- 7.898 1 +hi -tf +tc 



Fig. A41. — xcsao report of results of cross correla- 
tion of the object spectrum transforms against trans- 
forms of absorption and emission template spectra. 
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File: 1995 . 1122 . 0113 . 311 . 062165 .ms . imh HJD : 2450044.94573 



1995-NOV-23 10:37:37.00 
Object BCV: 25.298 
lambda: 3679.9 7551.7 
nbins: 4096x2 +tc 
Filter: 5 20 125 250 
Frac. endmask: 0.05 



Wavelength in Angstroms 
1 Corr . Template : FAST Emission Line Template 




Velocity [km/s] 

0.50-ht. peak fit, 9 pts . 



File: 1995. 1122. 0113. 311. 062165. ms. imh HJD : 2450044.94573 
Object:311. 062165 RA: 08:13:29.34 DEC : 29:02:12.00 2000.0 



1995-Nov-23 10:37:37.00 
Object BCV: 25.298 
lambda: 3800.0 6700.0 
nbins: 4096x2 +tc 
Filter: 5 20 125 250 
Frac. endmask: 0.05 



Wavelength in Angstroms 
2 Corr. Template: FAST Absorption Line Template 



- 1 1 1 


A 


- 1 1 1 


1114— 



Velocity [km/s] 



.50-ht. peak fit, 9 pts 



Fig. A44. — xcsao summary display for absorption 
cross-correlation. 



Fig. A42. — xcsao summary display for emission line 
cross-correlation. 



1995. 1122. 0113. 311. 062165. ms. imh JulDate: 2450044.94279 



1995-Nov-23 10:37:37.00 




5000 6000 7000 

Wavelength in angstroms 



Fig. A43. — xcsao summary display with labeled 
lines for emission line cross-correlation. 



a) VEL0CITY= 


14395.9542084403 


CZERR 


16.9255212706257 


VELQUAL = 




BCV 


25.2977199554443 


CZXC 


14395.9542084403 


CZXCERR = 


7.84048916091347 


CZXCR 


16.4074945116922 


BESTTEMP= 


' femtemp97' 


CZPARAMS= 


'2 4096 1 5 20 125 25 


XCTEMP1 = 


'femtemp97 14395.954 


XCTEMP2 = 


'fabtemp97 14414.752 


XCSAO 


'rvsao. xcsao 2.0b 05 


b) APVEL1 = 


' 14395 . 954 16 . 926 _ 2 


APVXC1 = 


'14395.954 7.840 16.4 


XT1_1 


'femtemp97 14395.954 


XT1_2 


'fabtemp97 14414.752 


XCSAO 


'rvsao. xcsao 2.0b 05 



7.840 16.4075 



182.443 0.000 
79 0.579 400.654 0.000 
:46 CZXC = 14395.95 R = 



000 C 
1.000 
16.41 



39:50 CZXC [1] ' 



Fig. A45. — Object spectrum header parameters 
added by xcsao: For a single spectrum, as shown in 
a), separate keywords are used for the results of the 
correlation. CZPARAMS gives the total number of 
templates, the number of points in the Fourier trans- 
forms, 2 if zero-padded, otherwise 1, and the four 
values for the cosine-bell filter used on the Fourier 
tranforms before they are correlated. For multi- 
spec spectra, as shown in b) VELOCITY, CZERR, 
VELQUAL, and BCV are combined on one line, as 
APVELn, and the best correlation results (CZXC, 
CZXCERR, CZXCR), best template spectrum se- 
quence number and name (BESTTEMP), and CZ- 
PARAMS are combined as APVXCn. In either case 
there is one line of information per template, giving 
template name, redshift velocity, velocity error, R- 
value, peak height (0.0 to 1.0), peak width at half- 
max (km/sec), template velocity and barycentric cor- 
rection, line-removal flag, and transform filter flag. 
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but options may be set to allow manual improvement 
of the line identifications and resulting redshift. The 
graphic cursor may be u sed to change fit and display 
parameters. Figure A46 shows the full parameter list 
for emsao. 

For each object spectrum file from the input list of 
spectra, and/or each aperture specified in the aper- 
ture list, the fitting subroutine, EMFIT, is called for 
the specified spectrum image band. Spectrum files 
are all read from the specified directory unless a full 
pathname is given in the spectrum list. Relative path- 
names may be used for spectra. If a directory is not 
set, the spectra are expected to reside in the current 
working directory. 

After the spectrum is loaded, it is renormalized, if 
the renormalize flag is yes. This should usually be 
done if the spectrum is in flux units. If the fixbad flag 
is yes, regions specified in the file named by badlines 
are replaced by straight lines. The spectrum is then 
smoothed nsmooth times. If objjplot is yes, the spec- 
trum is plotted and the plot, as shown in figure A47 is 
kept on the screen available for zooming and editing 
until a "q" is typed. 

To use a sky spectrum to compute noise statistics 
and to improve the error computations, an aperture 
or band must be specified. If skynum is not zero, a 
sky spectrum is read from that multispec aperture 
in the same file as specnum. If skyband is not zero, 
a sky spectrum is read from that multispec band in 
the same file. The sky spectrum, used to get the 
noise for error computations, is plotted if objjplot is 
yes. Previous results which have been saved in the 
spectrum image header may be used by setting linefit 
to no; in that case, all of the fitting below is skipped 
and the results are displayed. 

The continuum, computed by the IRAF curve fit- 
ting subroutine which is driven by the parameters set 
in the contpars pset task, is subtracted from the 
spectrum. If contsub-plot is yes, the spectrum is plot- 



ted with the continuum removed as in figure A48. 

The source of the initial redshift is specified by 
veLinit. If it is "guess" , the starting redshift is read 
from the parameter czguess. If it is "search", one 
line in the spectrum is identified by the program us- 
ing the table spe cified by the emsearch parameter, as 
shown in Figure |A49| . This table lists line centers in 
Angstroms and the wavelength range over which each 
one should be the strongest line. It can be modified 
by the user to match the data. The brightest line in 



skyband: 



esmooth 



logfile 
hardcop 



badlines . dat 



List of input spectra 

Spectrum aperture range if multispec or 2-D file 
Spectrum band if multispec 

Sky aperture range if multispec or 2-D file 
Sky band if multispec file 
Directory for input spectra 



fit 



bad li 



portions of spectrum 
Lines to eliminate from spectrum 
Renormalize spectrum before fitting (yes or no) 
Starting wave length of spectrum 
Ending wave length of spectrum 

Number of times to smooth spectrum for search and plot 
Initial velocity (search guess corr em file) 
Initial guess for Cz or z 

Specific cross-correlation template for initial velocity 

Wavelength to search around redshifted line center in angst 

Number of standard deviations above continuum to be line 

List of strongest emission lines for search 

Emission line list 

Directory for line lists 

Number of pixels to fit line peak 

Number of coeef f icients in line continuum fit 

Number of times to smooth spectrum for emission line fit 

If continuum is greater than this compute equivalent width 
Minimum fraction of mean line width for individual line 
Maximum fraction of mean line width for individual line 
Minimum equivalent width in sigma for individual line 
Velocity error if single line found (0 for gaussian fit err 
RMS dispersion error in angstroms 
Spectrum velocity correction 
Report mode (l=normal 2=one-line) 

Save results in binary archive record (yes or no) 
Save results in IRAF data file header (yes or no) 
Displays results on terminal 
List of log files 
Display output device 

Automatic hardcopy of results (yes or no) 
Plot results on terminal (yes or no) 
Hardcopy output device 

Graphical display mode (2=with line list 3=full screen) 
Velocity to plot 

Wait for cursor commands after plotting (yes or no) 

Label emission lines on spectrum plot (yes or no) 

Label absorption lines on spectrum plot (yes or no) 

Absorption line list 

Plot the object spectrum 

Plot the continuum-subtracted data 

Displays lots of (i.e. too many) intermediate results 
Graphics cursor input 



Fig. A46. — Parameters for emsao. 




5000 6000 
Wavelength in Angstroms 



Fig. A47. — emsao spectrum display if obj_plot is 
yes. 
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each region is assumed to be the one in the table, and 
its observed wavelength is saved. The redshift of the 
brightest of those lines is returned as an initial value 
to be refined by looking at more lines. 

If veLinit is "combination" , the initial redshift ve- 
locity is read from the spectrum header parameter 
VELOCITY; if "correlation" , it is read from the spec- 
trum header parameter CZXC; and if "emission", 
from CZEM. 

Search regions are read for each line in the file spec- 
ified by the emlines parameter. Figure |A50| shows the 
contents of such a file. Each region is then shifted by 
the guessed redshift and expanded in each direction 
by wspan Angstroms. All emission lines within each 
specified wavelength region are found. A spectrum 
pixel is assumed to be a line center if the pixel value 
is the max of the npfit neighbors on either side, and 
greater than linesig times the square root of the aver- 
age counts in those pixels. A second order fit is then 
made to the (2* npfit) +1 points centered on the peak 
to refine the center and peak height. The brightest 
line in each region is kept unless it has already been 
identified. Order matters the brightest line in a re- 
gion should be listed first, so that if it is the only one 
present in overlapping regions, it is correctly named. 

Before line profiles are fit, a copy of the spectrum is 
smoothed esmooth times using the same smoothing al- 
gorithm as is used before the line search is conducted. 
The parameter esmooth should be left at zero unless 
the data are especially noisy. It is best never to go 
above 2. Because local continuum values are impor- 
tant to the line fit, the continuum is removed from 
this copy of the spectrum using the same parameters 
as were used before conducting the line search. If 
contsub-plot is yes, the spectrum with the continuum 



removed is plotted as shown in figure A51 



Each identified line is checked to see if it is part 
of one of the groups of close emission lines listed in 
the file specified by emcombine, such as that shown in 
figure |A52| . If it is, all lines in the combination will be 
simultaneously fit. Those members of a line combi- 
nation which are not found are initialized at the red- 
shift of the most recently found line of the group. The 
missing line heights are assumed proportional to that 
line according to the relative heights in the emcom- 
bine file. The lines are fit by one to three Gaussians, 
along with an optional local continuum level using 
3 additional coefficients as set by nlcont. Redshift 
(computed from the wavelength of the center pixel co- 
ordinate), width, height, and errors are returned for 




5000 6000 
Wavelength in Angstroms 



Fig. A48. — The smoothed spectrum with the contin- 
uum subtracted is displayed by emsao if contsub_plot 
is yes. 



3727.30 3700.0 4000.0 Oil # 

4861.33 4840.0 4990.0 Hb # 

5006.84 4990.0 5400.0 Oil # 

6562.82 6500.0 7700.0 Ha # 



-2,100 to +22,000 km/sec 
-1,300 to +7,900 km/sec 
-1,000 to +23,500 km/sec 
-2,800 to +33,675 km/sec 



Fig. A49. — Emission lines for emsao's initial red- 
shift guess from the file emsearch.dat. The first col- 
umn is the line wavelength in Angstroms. The second 
and third columns are the minimum and maximum 
observed wavelengths to which it could be shifted. 
The fourth column is the name of the line. The com- 
ment field gives the redshift range over which the 
emission line will be in the given wavelength range. 
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Fig. A50. — Emission lines for emsao to look for 
in galaxy spectra from the file emlines.dat. If Z=0, 
OI at 6300 angstroms is from the earth's atmosphere. 
The first column is the line wavelength in Angstroms. 
The second and third columns are the minimum and 
maximum rest wavelengths from which to take con- 
tinuum values while looking for the line. The fourth 
column is the wavelength to fit on either side of the 
peak (for single line fits), and the last column is the 
name of the line. 




5000 6000 
Wavelength in Angstroms 



Fig. A51. — Unsmoothed spectrum with continuum 
removed, displayed by emsao if it contsub_plot is yes. 



each line. 

Each emission line is checked to see whether: 1) 
it has a fit center greater than zero, 2) it is wider 
than lwmin times the mean width of all of the iden- 
tified lines which are found and 3) narrower than Iw- 
max times that mean width, 4) it has an equivalent 
width more than Ismin times the error in the equiv- 
alent width, 5) it is too close to the blue edge of the 
spectrum, 6) it is too close to the red edge of the spec- 
trum, or 7) the error in the center of the Gaussian is 
zero. If dispmode is 1 or report-mode is 1, wavelengths 
and redshifts are printed for each line, with an X fol- 
lowed by the code for the test which failed at the end 
of entries which were omitted from the fit. If a line 
has been successfully fit, the rejection can be overrid- 
den interactively (if curmode is yes) using the + and - 
commands in cursor mode from the final display. If a 
line has been added or subtracted in cursor mode, a + 
or - at the end of the entry indicates that fact. A mean 
velocity is computed, weighted by the square of the 
error in the line centers and returned. If only a single 
line is found, the error is set to sigline, which should 
be set to the uncertainty in the dispersion function in 
Angstroms. This is usually significantly greater that 
the error in the fit to the center of the Gaussian. 

After all of the lines have been fit and a combined 
velocity has been computed, a zero-point redshift is 
computed by adding the solar system barycentric ve- 
locity correction, from a source specified by sveLcorr. 

The results are displayed as text to the devices 
specified by logfiles in the format specified by re- 
port-mode. Options include 1) one line per found 
emission line under a self-documenting summary, 2) a 
single line report listing the names of the lines which 
are found, and 3) a single line report listin g a veloc- 
ity for each searched-for line. Figure A53 shows the 
mode 1 report for the spectrum we are following. 



2 68 4958.91 

3 78 6548.06 
2 68 6716.440 



1.0 5006.84 1.0 0.00 0.0 # OI 

0.1 6562.82 1.0 6583.57 0.5 # N-Ha-N 

1.0 6730.815 -1.0 0.00 0.0 # S-S 



Fig. A52. — Emission line combinations for emsao 
to use with galaxy spectra from the file emcomb.dat. 
The first column is the number of lines to be fit to- 
gether. The second column is the additional wave- 
length range to fit blueward of the bluest peak and 
redward of the reddest peak. Then there are center 
wavelength and relative height pairs for each of the 
two or three lines to be fit together. 
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If displot is yes, the spectrum is plotted to de vice in 
the format specified by di spmo de. Figure A54 shows 
the mode 1 report; figure |A55| shows the mode 2 re- 
port. If hardcopy is yes, the same graph is automati- 
cally sent to a printer as well. 

If nsmooth is greater than zero, the displayed spec- 
trum is smoothed by a 1-2-1 sliding filter that many 
times. Absorption lines listed in the file ablines are 
labeled if dispabs is yes, and emission lines listed in 
the file emlines are labeled if dispem is yes. Both files 
are found in the directory linedir. 

If curmode is yes, the user can interact with the dis- 
play using the terminal cursor to zoom in on portions 
of the spectrum, identify lines and refit the emission 
lines, change the display format, edit the spectrum, 
or several other functions. 

If savejuel is yes, emission line redshift results are 
written into the spectrum image header in a form ap- 
propriate to the spectrum format: two entries plus 
one per line if multispec; otherwise one, value per 



keyword, as shown in Figure A56 



A. 3. Creating a Spectrum with linespec 

linespec is an IRAF task for making a spectrum 
from a list of emission and/or absorption lines which 
is driven by the parameters shown in Figure A57. It 
was written to create templates for use by xcsao. An 
example of such a line list, with center wavelengths, 
widths, and heights for each line whic h may appear in 
the spectrum, is shown in figure |A58| . The same line- 
defining parameters are used as are output of emsao, 
so low-noise templates can be easily made. The line 
list is read from the file specified by the linefile pa- 
rameter in the directory designated by the parameter, 
linedir. If linedir is null, the file is assumed to be in 
the current working directory. 

A blank (all zero) spectrum, with the object name 
specobj and file name specfile in the directory specdir 
is created. It is linear in wavelength, with a resolution 
in Angstroms given by pixJambda, and a range from 
stJambda to end-lambda. Spectral world coordinate 
system information is written to the header, with the 
standard FITS parameters CRPIX1, CRVAL1, and 
CDELT1. If verbose is yes, files specified by logfiles 
are opened and a header is written. 

The center of each line in the table is redshifted 
according to the parameters zspec and velspec. vel- 
spec (cz=apparent Doppler shifting velocity) is used 
unless zspec (AA/A) is not zero, in which case zspec is 



r2> emsao 1995 . 1122 . 0113 . 311 . 062165 .ms 

IRAF rvsao. emsao 2.0b minkgcfaO 06-Oct-1997 16:47 

File: 1995 . 1122 . 0113 . 311 . 062165 .ms [1] 

Object: 311.062165 RA: 08:13:29.34 Dec: 29:02:12.0 2000.0 
Observed 23-Nov-1995 10:37:37.00 = JD 2450044.9428 BCV: 0.00 
Combined vel = 14365.86 +- 19.93 km/ sec, z= 0.0479 
Emission vel = 14365.86 +- 13.12 km/ sec, z= 0.0479 for 9/10 lines 



Line 
Oil 



OIII 
OIII 



Rest lam 
3727.30 
4861.33 
4958.91 
5006.84 
6300.23 
6548.06 
6562.82 
6583.57 
6716.44 



Obs . lam 
3905.71 
5094.26 
5196.17 
5247.03 
6602.31 
6861.83 
6877.34 
6898.86 
7038.27 



Pixel 
154.60 
963.18 
1032.51 
1067.11 
0.00 
2165.68 
2176.23 
2190.87 
2285.71 



0479 
0479 
0478 
0480 
0000 
0479 
0479 
0479 
0479 



vel 
14350.03 
14364.32 
14343.84 
14381.68 
0.00 
14365.28 
14367.50 
14357.17 
14365.01 



dvel 
76.84 

8.95 
88.51 
36.56 

0.00 

9.32 



eqw 
7.12 
7 .50 
0.53 
2.40 
0.00 
6.59 
33.81 
17.35 
6.00 



0.002 

0.000 XI 

0.029 

0.709 

0.177 

0.031 



Fig. A53. — emsao report for an emission line galaxy 
when report-mode is 1. For each emission line, the 
name, rest wavelength, observed wavelength, pixel 
number at the observed wavelength, redshift in frac- 
tion of the speed of light and km/sec and error in 
km/sec, equivalent width, and relative weighting fac- 
tor are printed. An "X" , followed by a rejection code, 
after the last column indicates that the line was not 
used. 



395-NOV-23 10:37:37. OC 



Object:311. 062165 RA: 08:13:29.34 DEC: 29:02:12.00 2000.0 

I ' ' ' ' I ' ' ' ' I ' ' ' 'I 




n r F r i i *f i i i i 



5000 6000 7000 

Wavelength in angstroms 



Fig. A54. — The emsao summary display shows the 
spectrum and line fit results if dispmode is 1. 
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1995. 1122. 0113. 311. 062165. ms 



vel= 14365.86 (13.13) km/sec 








i a 



5000 6000 
Wavelength in angstroms 



Fig. A55. — The emsao summary display shows only 
the spectrum if dispmode is 2. 



VELOCITY= 
CZERR 
VELQUAL = 
BCV 
CZEM 

CZEMERR = 
CZEMNL = 
CZEMNLF = 



EMLINE4 = 



EMLINE9 = 
EMLINE10- 



14365.8563351488 
19.9269054499303 



0. 



14365.8563351488 
13.1179861568178 
10 
9 

'Oil 3727.30 3905.71 136.27 2.5747 14350.03 76.84 7.12 0.000' 
'Hb 4861.33 5094.26 742.95 1.9504 14364.32 8.95 7.50 0.031' 
'OIII 4958.91 5196.17 65.501 1.7572 14343.84 88.51 0.53 0.000' 
'OIII 5006.84 5247.03 191.85 2.6152 14381.68 36.56 2.40 0.002' 
'OI 6300.23 6602.31 0.0000 0.0000 0.00 0.00 0.00 0.000' 
'Nl 6548.06 6861.83 774.78 2.3900 14365.28 9.32 6.59 0.029' 
'Ha 6562.82 6877.34 3955.1 2.3902 14367.50 1.87 33.81 0.709' 
'N2 6583.57 6898.86 1956.0 2.4612 14357.17 3.75 17.35 0.177' 
'SI 6716.44 7038.27 675.48 2.4088 14365.01 8.90 6.00 0.031' 
'S2 6730.81 7053.84 535.17 2.3420 14387.80 11.05 4.60 0.020' 
'rvsao. emsao 2.0b 06-Oct-1997 17:25 CZEM = 14365.86 LINES = 9' 



Fig. A56. — These keywords save the emsao re- 
sults in the spectrum's FITS header. For each emis- 
sion line, the name, rest wavelength, observed wave- 
length, observed height and width, redshift and error 
in km/sec, equivalent width, and relative weighting 
factor are saved. 
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Fig. A58. — List of emission lines from linespec 
input file. The first column is the center wave- 
length in Angstroms, the second is the half- width in 
Angstroms, the third is the height in arbitrary units 
to be used in the output spectrum. The last column 
identifies the line. 



linefile= 
(linedir= 
(linewid= 
(maxwidt= 
(zspec = 



(end lam= 

(spec_pl= 
(spec_in= 
(device = 
(plotter= 
(verbose= 
(logfile= 



n or absorption lines 

2.) HWHM of instrument resolution in Angstroms 

no) Use maximum of line or instrument width or convolve (yes or no) 

0.) Shift line centers to this dlambda/lambda (overrides velspec) 

250000.) Shift line centers to this velocity 

0. ) Continuum level for created spectrum 
QSO template) Object name for created spectrum 

qsotemp Filename for created spectrum 

) Directory for created spectrum 

2000.) Starting wave length to save 

10000.) Ending wave length to save 

1 . ) Wavelength per pixel in Angstroms 
yes) True to plot the created spectrum 
yes) True to interact with the created 

stdgraph) Display device 

stdplot) Hardcopy output device 

yes) Print summary to log file (yes or not) 

STDOUT, linespec. log) List of log files 

no) Displays intermediate results 

) Graphics cursor input 



Fig. A57. — Parameter list for linespec 
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used for the redshift. The linewidth, if it is tabulated 
in kilometers per second, is converted to Angstroms 
at the shifted line center. The line width is also 
broadened appropriately if the line is redshifted. The 
redshift velocity is put in the spectrum header us- 
ing the VELOCITY keyword. If maxwidth is yes, the 
linewidth parameter is used as the width of the line if 
it is greater than the tabulated width. If maxwidth is 
no, the width from the table is used, and the linewidth 
smoothing is done later. For each line, a Gaussian at 
the shifted center wavelength, half-width, and tabu- 
lated height is added to the spectrum. As each line is 
added to the spectrum, the line definition is written 
to the spectrum's he ader, producing a table such as 
that shown in figure A59| . After all of the lines are 
computed, a continuum level specified by the param- 
eter continuum is added to the spectrum. 

The spectrum is p lotted to a graphics device, as 
shown in figure |A60| , if spec-plot is yes. After the 
graph is displayed to the terminal, if specAnt is yes, 
the user must interact with it using single character 
cursor commands, such as the "z" for zoo m co mmand, 
the result of which is shown in figure A61 . If the 
"@" command is given in this mode, the graph is 
immediately sent to the device specified by plotter. 
The "q" command must be typed to leave the graph 
and proceed. 

If maxwidth is no, the resolution of the spectro- 
graph is simulated by convolving the entire spectrum 
with a Gaussian of height 1.0 and sigma (half width 
at half of maximum) of linewidth. If specjplot is yes, 
the spectrum is displayed again. 

The linespec version and the date the program is 
being run are written to the spectrum image header, 
the spectrum is written, and the file is closed. 



EMLINEl = 
EMLINE2 = 
EMLINE3 = 
EMLINE4 = 
EMLINE5 = 
EMLINE6 = 
EMLINE7 = 
EMLINE8 = 
EMLINE9 = 
EMLINEl 0= 
EMLINEl 1= 



Mgll 
Oil 



2798 . 
3727 . 
4861 . 
4958 . 
5006. 
6300 . 
6548 . 
6562 . 
6583. 
6717 . 
6731 . 



00 2798.00 20.00 2. 

30 3727.30 60.00 2. 

33 4861.33 115.00 2. 

91 4958.91 50.00 2. 

84 5006.84 127.00 2. 

23 6300.23 32.00 2. 

06 6548.06 67.00 2. 

82 6562.82 700.00 2. 

57 6583.57 194.00 2. 

00 6717.00 132.00 2. 

30 6731.30 100.00 2. 
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90 
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70 
40 



60 
60 



145.38 
436.15 
720 . 66 
325.86 
859.52 
192.51 
453 .45 
4562 .06 
1264.34 
893.36 
676.79 



, 000' 
. 000' 
, 000' 
. 000' 
, 000' 
. 000' 
, 000' 
. 000' 
. 000' 
, 000' 
. 000' 




5000 6000 
Wavelength 



Fig. A60. — Graph of an emission line spectrum pro- 
duced by linespec. 




Wavelength 



Fig. A59. — These multiple-valued keywords are 
added to the header of the synthesized spectrum 
to show what lines were added. For each emission 
line, the name, center, redshifted center, height, half- 
width, velocity, and equivalent width. 



Fig. A61. — Graph of hydrogen alpha region of the 
same spectrum displayed using the "z" cursor com- 
mand in linespec. 
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A. 4. Creating a Composite Spectrum with 
sumspec 

sumspec is an IRAF task for making a compos- 
ite spectrum which is adjusted for velocity. It re- 
bins input spectra to a common wavelength or log- 
wavelength range and resolution, with optional con- 
tinuum removal and normalization. While it was orig- 
inally intended for making templates to be used by the 
xcsao cross-correlation task, it has turned out to be 
very useful for rebinning individual spectra as well as 
combining multiple spectra. The parameter list for 
sumspec is shown in Figure A62 . 

The spectra to be added are specified by the list 
of files in the parameter spectra and/or the list of 
spectrum numbers in specnum. Unless a pathname is 
specified as part of the filename in spectra, each spec- 
trum is expected to be found in the directory specdir. 
All spectra are taken from the band, specband, in mul- 
tispec spectra. Files specified by logfiles are opened to 
receive logging information. If debug is yes, additional 
information about the input files and the progress of 
the program is written to STDERR. 

If either stJambda or endJambda is set to INDEF, 
all of the spectrum headers are read to find the lim- 
its of the overlap of their wavelength coverage, and 
the missing limit (s) of the output spectrum are set 
accordingly. 

Before the first spectrum in the list is added, a 
blank (all zero) spectrum, with the object name com- 
pobj and file name compfile in the directory com- 
pdir is created. If complog is yes, it is linear in log- 
wavelength; otherwise, it is linear in wavelength. Ei- 
ther way, the output spectrum contains npts data 



( specnum= 
compf ile= 



(normin = 
(fixbad = 
(badline= 
(linedir= 



(cont_pl= 

(end_lam= 

(npts 

(complog= 



(spec_pl= 



(plotter 
(logfile 



List of input spectra 

Spectrum aperture range if multispec or 2-D file 

Spectrum band if multispec file 

Directory for input spectra 

Filename for output composite spectrum 

Object name for output composite spectrum 

Directory for output composite spectrum 

Save names of input files in output file header 

Renormalize spectrum before processing (yes or no) 

Eliminate portions of spectrum with bad lines (yes or no) 

Lines to eliminate from spectrum 

Directory for line list 

Remove continuum from input spectra (subtract , divide, no) 

Remove lines from input spectra (yes, no) 

Input absorption line rejection in sigma of fit 

Input emission line rejection in sigma of fit 

Plots the continuum- subtracted data 

Continuum level added to composite spectrum 

Starting wave length to save 

Ending wave length to save 

Number of pixels for output spectrum 

Rebin into log wavelength (yes or no) 

Rebin interpolation mode 

Renormalize spectrum after rebinning (yes or no) 
True to plot each spectrum 

True to interact with each spectrum graph 
True to plot the composite spectrum 
True to interact with the composite spectrum 
Redshift composite velocity in km/sec 

Redshift composite velocity in dlambda/lambda (overrides 

Spectrum velocity correction (none, file , helio, bary, hfile) 

Smooth displayed spectrum (0=no else number of times) 

Display device 

Hardcopy output device 

List of log files 

Displays intermediate results 

Graphics cursor input 



points between stJambda, or its log, and endJambda, 
or its log. World coordinate system information is 
then written to the header. 

Spectra are read from the input list one at a time 
and renormalized if renormalize is set to yes. This 
should be done if the spectra are fluxed or whenever 
pixel values are much bigger or smaller than their 
variation. If specjplot is yes, each spectrum is plot- 
ted to the devic e specified by the parameter device, 
as in figure |A64 If specJnt is yes, the user can in- 
teract with that graph using single character cursor 
commands and must type a "q" to proceed with the 
task. If the "@" command is given in this mode, the 
graph is sent to the device specified by plotter. 

Each spectrum is rebinned using an interpola- 
tion mode specified by interpjmode into a npte-pixel 
wavelength-linear spectrum covering the range com- 
puted above. If complog is yes, the rebinned spectrum 
is linear in log- wavelength; if complog is no, the re- 
binned spectrum is linear in wavelength. If specplot 
is yes, the rebinned spectrum also is plotted. 

If s-contin is set to subtract or divide instead of no, 
the IRAF curve fitting subroutines, are used to fit a 
continuum to each input spectrum. The CONTSUM 
task, described in the following section, sets all of the 
appropriate parameters. 

As it is rebinned, each input spectrum is redshifted 
according to the parameters ztemp (z — A A/ A), if it is 
not zero, or veltemp (cz = apparent Doppler shifting 
velocity). A velocity correction to the solar system 
barycenter is removed according to the sveLcorr pa- 
rameter. Set it to none if the input spectra have not 
been shifted. If "file", BCV is used if present in the 
file header, or else HCV. If "hfile", the header pa- 
rameter HCV is always used. If neither is found, no 
correction is made. If "heliocentric" or "barycentric" 
corrections are chosen, position and time parameters 
are read from the spectrum data file header, and the 
correction is computed as described in section |A.6| be- 
low. 

Emission and/or absorption lines may be removed 
from each input spectrum if reject is yes. If the 
spectrum header parameter SUBCONT is present, its 
value overrides that of reject. The parameters absrej 
and emrej set the lower and upper acceptable limits 
for input spectrum pixels in standard deviations of 
the continuum fit to the spectrum. 

Graphs of the continuum-removed, apodized data 



Fig. A62. — Parameter list for sumspec. 



are displayed, as in figure A64 if the contjplot param 
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eter is yes. 

The composite spectrum is plotted to the d evice 
specified by the parameter device, as in figure |A65| , 
if compjplot is yes. If comprint is yes, the display is 
held to allow the user to interact with it using single 
character cursor commands. If the "@" command is 
given in this mode, the graph is sent to the device 
specified by plotter. A "q" command allows the task 
to proceed. 

After each spectrum is added, the sumspec ver- 
sion and the current date are written to the output 
spectrum image header and that spectrum image file 
is written. After the last spectrum is added the file 
is closed. Figure |A66| shows the parameters added 
to the spectrum header. EXPTIME is the total ex- 
posure time from all of the input images. A one- 
dimensional image is produced, so DISPAXIS is al- 
ways 1. DC-FLAG is 1 if the output spectrum is 
log- wavelength, in which case the log- wavelength of 
the first pixel is given in both CRVAL1 and WO, 
and the log- wavelength per pixel is given in CDELT1 
and WPC. The VELOCITY is set by the velcomp or 
zcomp parameter. If no velocity is specified, it is left 
at 0. If savenames is yes, the instrument, filename, 
input redshift velocity and barycentric velocity cor- 
rection are written to the header. 

A. 5. Removing a Spectrum's Continuum us- 
ing contpars or CONTSUM 

contpars and CONTSUM are IRAF parameter 
files which are used by a single subroutine to fit and 
remove a continuum from object or template spectra. 
The subroutine uses the IRAF interactive curve fit- 
ting subroutine package and is based on the technique 
used in the IRAF RV package (Fitzpatrick 1993). 
The parameters which are set are shown in figure 
A67. The same subroutine is called prior to cross- 



correlation by xcsao, an emission line search by em- 
sao, and before or after summation of spectra by 
sumspec. If the SUBCONT header keyword in a 
template spectrum is set to F, no continuum is fit or 
removed from that template spectrum. If the keyword 
is T or not present, the continuum is fit and removed. 

For a spectrum, such as that in figure |A68| , the 
continuum is fit using the IRAF curve fitting sub- 
routine. It may be used interactively by setting the 
C-interactive parameter to yes. In that case, a graph 
of the fit an d t he sp ectrum is displayed as shown in 
figures |A69| and |A70| . The result may be examined or 



a) 




b) 




Fig. A63. — Input spectra for sumspec. 



b) 4 




Fig. A64. — Rebinned sumspec input spectra with 
the continuum removed. 
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Fig. A65. — Output composite of two input spectra 
plotted by sumspec. 



EXPTIME = 

DISPAXIS= 

DC-FLAG = 

CRPIX1 = 

CRVAL1 = 

CDELT1 = 

WO 

WPC 

NP1 

NP2 

VELOCITY= 

HISTORY 

HISTORY 

HISTORY 

HISTORY 



1 
1 . 

3.55497165126008 
1.47901535349635E-4 

3.55497165126008 
1 . 47901535349635E-4 

1 

2048 
0. 

RVSAO. SUMSPEC 2.0 run 05-May-1997 15:14 
Add FAST 1 9 95 . 4 02 . 1 1 8 . 1 32 65p0 12 . ms at CZ 6446.84 - 
Add FAST 1 9 95 . 1 1 22 . 1 1 3 . 3 1 1 . 62 1 65 . ms at CZ 14385.28 
Add FAST 1 9 9 6 . 1 23 . 052 . 4 65 . 1 1 857 . ms at CZ 11585.76 



1.07 
25.30 
-28 .74 



Fig. A66. — These parameters are set in the output 
spectrum produced by sumspec. 



refit after individual pixels are removed. A "q" com- 
mand is needed to exit from the interactive graph. 

The cjunction parameter specifies the function to 
be used to fit the continuum. If the continuum is 
irregular, it is best to use the spline3, but splinel, 
Legendre, or Chebyshev polynomials may be used as 
well. The parameter order sets the order of the fit. 
The portion of the the spectrum to be fit is specified, 
in image section format, by the parameter sample, 
Before the fit, the spectrum is averaged in groups of 
naverage pixels. 

Pixels with values outside of the limits sdow -reject 
and sJiigh-r eject, for spectra, and tdowjreject and 
t-high-reject, for templates are rejected. Separate lim- 
its are used because composite template spectra will 
have much higher signal-to-noise levels than object 
spectra and limits may need to be set differently to 
fit a good continuum. These limits are specified in 
standard deviations of the fit. To make sure that the 
wings of lines are thoroughly eliminated, this rejec- 
tion process is repeated niterate times. Each time, 
a number of pixels, specified by the grow parameter, 
are rejected on each side of the rejected pixel. 

Normally, the continuum which is fit is subtracted 
from the spectrum, preserving the signal-to-noise of a 
spectrum when it is dominated by photon noise. To 
remove the continuum of both template and object 
spectra by division, the template spectrum header 
keyword DIVCONT is set to T. In that case, noise 
in portions of the spectrum with low signal-to-noise 
levels will be amplified. If ems ao.cont sub -plot, xc- 
sao.contsubjplot, or sumspec. contjplot is yes, the spec- 
trum is plotted with the continuum removed, as 



shown in figure A71 



A.6. 



Correcting Radial Velocities to the So- 
lar System Barycenter 



_high_= 
~high_= 



flpar = 



*) Sample of points to use in fit 

ie3) Fitting function 

(spline3 [ legendre I chebyshev | splinel ) 

1) Number of points in sample averaging 

10) Order of fitting function 

2.) Spectrum continuum low rejection in sigma of fit 

2.) Spectrum continuum high rejection in sigma of fit 

2.) Template continuum low rejection in sigma of fit 

2.) Template continuum high rejection in sigma of fit 

10) Number of rejection iterations 

no) flush pfile on assign? 

qD 

0) 



Fig. A67. — contpars parameters as shown by IRAF 
eparam task. 



To compare redshift velocities observed when the 
earth is at different positions in its orbit, the velocity 
of the sun relative to the earth is added to the red- 
shift velocity. This heliocentric velocity correction is 
more accurately done to the true stationary point of 
the solar system barycenter, its center of mass. This 



barycentric velocity correction is shown in figure A72 
A second, smaller correction is added for the motion 
of the observer relative to the center of the earth as 
it rotates. This vector velocity is projected in the di- 
rection of the observed object, and that component is 
saved. 
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NO AO/I RAF V2.11 EXPORT mink@xena Thu 16:11:15 08-Jan-98 
func=spline3, order=15, low_rej=2, high_rej=2, niterate=10, grow=1 
total=2635, sample=2635, rejected=1061, deleted=0, RMS= 22.31 



5000 6000 
Wavelength in Angstroms 




Fig. A68. — Input spectrum before the continuum is 
removed according to contpars. 



Fig. A70. — Here is a closeup of the region around 
Ha, showing the continuum which has been fit and 
giving a better view of the rejected pixels. 



NOAO/IRAF V2.11 EXPORT mink@xena Thu 16:10:19 08-Jan-98 
func=spline3, order=15, low_rej=2, high_rej=2, niterate=10, grow=1 
total=2635, sample=2635, rejected=1061, deleted=0, RMS= 22.31 




L 



5000 6000 
Wavelength in Angstroms 



Fig. A69. — On this display shown during interactive 
continuum fitting, spectrum pixels rejected from the 
fit are marked with diamonds. 



Fig. A71. — Figure |A68| spectrum with its continuum 
removed. 



39 




Fig. A72. — The correction of a radial velocity ob- 
served from the earth to that which would be observed 
from the center of mass of the solar system is shown 
at two different times in the earth's orbit. 

The xcsao, emsao, and sumspec tasks of the 
RVSAO package compute this correction using sub- 
routines which read the time of observation, object 
position, and observatory position from the spectrum 
image header. Several common alternative keywords 
are built into those subroutines. RA, DEC, and 
EPOCH give the right ascension, declination, and 
equinox of the observed object. SITELONG, SITE- 
LAT, and SITEELEV for longitude, latitude, and al- 
titude, respectively give the observer's location on the 
earth. OBS-DATE yields the observation date and 
UTMID, the midtime of the observation. If UTMID 
is not found in the header, UT, assumed to be the 
end time of the observation, or UTOPEN, the start 
time of the observation are used in combination with 
EXPOSURE or EXPTIME, which give the duration 
of the observation in seconds to compute a midtime. 
The requested heliocentric or barycentric correction 
is then calculated at that midtime. 

Since there is really no standard for the meaning of 
these keywords, a separate task, bcvcorr, has been 
added to RVSAO to allow several alternate ways of 
specifying these three major pieces of information, 
bcvcorr can write its result to the header of the im- 
age which it is processing; the other RVSAO tasks 
will use this value when their sveLcorr parameter is 



set to "file". First, the sky direction of the object is 
set. The right ascension, declination, and equinox of 
the object's sky position are read from the keywords 
specified by the keyra, keydec, and keyeqnx parame- 
ters. If those are not all found, the position is read 
directly from the parameters ra, dec, and equinox. 

The time for the velocity correction is read from 
a Julian date keyword specified by keyjd. If none 
is found, the date is taken from the header keyword 
specified by keydate. The UT midtime is taken from 
the keyword in key mid. If that is not found, the obser- 
vation midtime is computed by finding the exposure 
duration from the keyexp keyword value and adding 
half of it to the start time from the keystart keyword 
value or subtracting half of it from the end time in 
the keyend keyword value. If no time can be found 
in the header, the Julian date is read directly from 
the gjd parameter, if it is greater than zero, or the 
midtime is taken from the year, month, day, and ut 
parameters. 

The heliocentric Julian Date (HJD) is the time at 
which the light from the object for this observation 
reachs the sun. It is needed if multiple radial velocity 
observations of an object are to be compared accu- 
rately. The HJD is not used in computation of the 
velocity correction, which is dependent on the time 
of the observation at the earth. It may be read from 
the header keyword keyhjd or set by the parameter 
hjd (only if > 0). If none is set, the HJD is computed 
from the observation time. 

If the obsname parameter is "file" , the observatory 
name and position is read from the image header us- 
ing the keywords in keyobs, for the name, keylat, for 
the latitude, keylong, for the longitude, and keyalt, for 
the altitude. Otherwise, the value of obsname is used 
to get a position from IRAF's observatory database. 
If the string is not found there, the longitude, latitude, 
and altitude are given by the parameters obslong, ob- 
slat, and obsalt. 

If verbose is yes, the information in Figure |A73| is 
printed to standard output. If savebcv is yes, sev- 
eral spectrum image header keywords are set. BCV 
is the barycentric velocity correction in km/sec. The 
midtime of the observation is stored three ways: UT- 
MID, the Universal Time, GJDN, the geocentric Ju- 
lian Date, and HJDN, the heliocentric Julian Date. 
These times are displayed by xcsao and emsao as 
the time of observation. 
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Fig. A73. — bcvcorr parameters as shown by IRAF 
eparam task. 

RA: 8:13:29.3, Dec: 29:02:12.0 2000.0 

Using default observatory database: noao$lib/obsdb . dat 
Using observatory parameters for database entry: flwo 
flwo lat 31:40:51.4 , long 110:52:39.0, alt 2320.0 
Julian date is 2450044.94279 at 11/23/1995 10:37:37.0 UT 
Object at ra 8:13:29.340 dec 29:02:12.00 eq 2000 
Heliocentric Julian date: 2450044.94573 
gbcvel = 25.2207 ghcvel = 25.2165 geovel = 0.0770 
bcv = 25.2977 hcv = 25.2936 computed 
bcv = 0.0000 hcv = INDEF from file 



Fig. A74. — Output from bcvcorr 
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